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ABSTRACT 

A  quantum  mechanical  method  is  used  to  investigate  the 
linear  and  nonlinear  collisional  processes  in  an  unmagnetized 
plasma,  which  arise  out  of  the  scattering  of  two  particles 
through  an  effective  Coulomb  field  via  emission  or  absorption 
of  one  and  two  plasmons. 

The  validity  of  Kadomtsev's  assumption  (1965)  concern¬ 
ing  the  ion  nonlinear  Landau  damping  has  been  questioned 
(Sloan  and  Drummond,  1970) ,  inasmuch  as  the  electron  nonlin¬ 
ear  wave-particle  effects  make  the  wave  further  destabilizing 
The  nonlinear  dielectric  function  and  its  effect  on  wave- 
particle  coupling  constant  is  investigated  and  it  is  shown 
that  a  reduction  in  the  growth  rate  occurs  and  that  the  ampli 
tude  dependent  frequency  shift  causes  stabilization. 

The  four-wave  interactions  and  the  amplitude  dependent 
frequency  shift  for  electromagnetic  waves  (light,  whistler 
and  Alfven)  propagating  parallel  to  an  external  magnetic 
field  are  studied  by  employing  the  diagrammatic  method  of 
field  theory. 
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CHAPTER  I 

INTRODUCTION 


The  linear  theory  in  plasma  physics  had  been  adequate 
for  the  understanding  of  small  amplitude  plasma  oscillations 
and  indeed  it  was  believed  that  the  problem  of  the  stable 
confinement  of  plasma  could  be  solved  within  this  linear 
theory.  More  recent  theoretical  development  and  experimen¬ 
tal  measurements  have  demonstrated  a  whole  chain  of  new 
plasma  instabilities  and  in  fact  the  presence  of  instabili¬ 
ties  is  the  most  characteristic  attribute  of  this  state  of 
matter. 

The  fundamental  question  as  to  how  the  instabilities 
develop  and  whether  or  not  they  play  a  significant  role  in 
plasma  phenomena  can  only  be  answered  by  nonlinear  theory. 

In  the  past  few  years,  the  application  of  nonlinear  theory 
(Sturrock,  1957;  Ginzburg  and  Zhelemiakov,  1958;  Kadomtsev, 
1965  and  Vedenov,  1965)  has  undergone  such  vigorous  develop¬ 
ment  as  to  result  in  the  formulation  of  some  clear  physical 
concepts  regarding  these  nonlinear  mechanisms.  However 
while  the  achievements  of  the  above  authors  were  notable,  the 
range  of  nonlinear  interactions  was  by  no  means  completely 
understood  nor  was  a  systematic  methodology  developed.  In 
addition  to  the  class  of  problems  in  which  the  amplitude  of 
the  wave  fields  is  of  primary  concern,  nonlinear  theory  can 
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be  used  to  explain  mode  conversion.  This  process  is  of  funda¬ 
mental  importance  since  mode  conversion  leads  to  electromag¬ 
netic  radiation  by  which  interstellar  and  solar  plasmas,  for 
example,  can  be  studied. 

When  considering  the  problem  of  nonlinear  interactions, 
the  quantum  mechanical  calculations  are  more  straightforward 
and  less  difficult  than  the  corresponding  classical  calcula¬ 
tions.  One  of  the  most  successful  early  methods  of  attacking 
this  problem,  which  is  still  used  extensively  today,  is  the 
canonical  transformation  technique.  However,  one  would  prefer 
a  more  systematic  method  of  obtaining  system  properties  and 
it  was  shown  in  a  series  of  revolutionary  papers  (Vedenov, 
et  al. ,  1962,  1963;  Pines  and  Schrieffer,  1962)  that  quantum 
field  theory,  previously  restricted  to  elementary  particle 
physics,  provides  such  a  systematic  method.  An  essential 
role  in  the  field  theoretic  treatment  of  the  many  body  prob¬ 
lem  is  played  by  the  Green's  function  or  propagator  of  the 
system.  The  propagators  can  be  calculated  relatively  more 
easily  by  using  the  well  known  Feynman  diagram  technique 
which  has  the  additional  advantage  of  providing  a  pictorial 
representation  corresponding  to  each  term  in  the  perturbation 
expansion . 

From  the  above  point  of  view,  the  waves  in  a  plasma 
can  be  thought  as  being  composed  of  quasiparticles  (quanta  of 
waves) .  These  quasiparticles  interact  with  the  particles  of 
the  plasma  and  with  each  other  and  their  interactions  are 
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described  in  terms  of  an  interaction  Hamiltonian.  The  inter¬ 
action  Hamiltonians  for  particles  and  plasmons  (quanta  of 
plasma  oscillations)  and  particles  and  phonons  (quanta  of  ion 
sound  waves)  were  obtained  by  Pines  and  Schrieffer  (1962). 

They  used  these  together  with  the  "Fermi  Golden  Rule"  to 
write  the  quasilinear  equations.  Pines  and  Schrieffer, 

(1962)  and  Vedenov,  et  al.  (1962)  pointed  out  that  the  Landau 
damping  or  growth  process  which  forms  the  basis  of  quasilin¬ 
ear  theory,  can  be  described  as  the  competition  between  ab¬ 
sorption  and  stimulated  emission  of  quantized  plasma  waves 
or  plasmons.  Similarly,  the  quantum  form  of  the  equation 
describing  nonlinear  Landau  damping  is  discussed  by  Rosen- 
bluth,  et  al. ,  (1969);  Ross,  (1969);  Harris,  (1969)  and 

Krishan  and  Fukai,  (1971). 

The  rate  of  change  of  particle  distribution  is  de¬ 
rived  by  Wyld  and  Pines  (1962).  They  assumed  that  the  ma¬ 
trix  element  for  a  Coulomb  collision,  4iTe2/q2  must  be  modified 
by  the  factor  e”1(q,a))  (Thompson,  1962)  where  e  is  the  di¬ 
electric  function  of  the  plasma,  fiq  is  the  momentum  transfer, 
and  Kw  is  the  energy  transfer  in  the  collision.  The  equa¬ 
tion  of  Wyld  and  Pines  reduces  to  the  well  known  Boltzman 
equation  in  the  classical  limit.  Quantum  mechanical  calcu¬ 
lations  for  f our-plasmons  (Langmuir  waves)  in  an  unmagnetized 
plasma  have  been  carried  out  by  Zakharov  (1967)  and  for  three 
plasmon  interactions  for  longitudinal  and  transverse  waves  by 
Krishan  (1968)  and  by  Krishan  and  Selim  (1968) .  Attention 
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should  also  be  directed  to  the  paper  of  Gailitis,  et  al. 

(1966)  in  which,  although  the  calculations  are  classical, 
the  language  is  quantum  mechanical.  This  deals  with  the 
interaction  of  plasmons,  phonons  and  photons  in  an  isotropic 
plasma.  A  quantum  mechanical  theory  of  nonlinear  phenomena 
in  a  very  strong  magnetic  field  was  developed  by  Walters  and 
Harris  (1967,  1968).  The  three  plasmon  interaction  of  Wal¬ 
ters  and  Harris  had  previously  been  derived  classically  by 
Aamodt  and  Drummond  (1964). 

In  chapter  II,  the  occupation  number  formalism  (sec¬ 
ond  quantization)  and  the  motion  of  charged  particles  in  a 
uniform  magnetic  field  are  introduced.  The  electromagnetic 
field  in  a  plasma  is  quantized  as  has  previously  been  done  by 
Kihara,  Aono,  and  Dodo,  (1962);  Alekseev  and  Nikitin  (1966) 
and  Harris  (1969) .  The  wave-particle  and  particle-particle 
interaction  Hamiltonians  are  derived  and  represented  by 
first  order  Feynman  diagrams.  The  theory  of  the  scattering 
matrix  and  its  representation  by  Feynman  graphs  are  developed 
and  represented  in  a  suitable  form  for  studying  the  nonlinear 
plasma  interaction.  General  rules  for  the  calculations  of 
the  matrix  elements  (coupling  constants)  for  higher  order 
diagrams  are  established. 

Chapter  III,  deals  with  the  transition  matrix  elements 
for  the  scattering  of  two  particles  through  effective  Coulomb 
field  and  the  emission  (absorption)  of  one  and  two  plasmons 
(waves)  in  an  unmagnetized  plasma.  Using  these  matrix  elements 
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and  the  "Fermi  Golden  Rule"  one  obtains  the  linear  and  non¬ 
linear  collisional  damping  for  longitudinal  plasma  waves. 

In  chapter  IV,  the  dielectric  function  and  linear 
Landau  growth  rate  for  ion  acoustic  waves  are  derived  quan¬ 
tum  mechanically,  in  such  a  way  as  to  emphasize  the  similarity 
in  the  classical  and  quantum  mechanical  derivations.  Using 

the  same  method,  the  nonlinear  dielectric  function  e„_  is 

NL 

derived  and  its  effect  on  the  wave-particle  coupling  constant 
is  discussed.  Solving  e  =  0/  the  frequency  shift  and  its 
stabilization  influence  on  the  unstable  waves  is  discussed. 

Using  the  modified  coupling  constant  and  the  Fermi  Golden  Rule, 
the  quasilinear  equations  for  the  plasma  particles  are  obtained. 
These  equations  yield  the  ratio  of  the  heating  rates  for  the 
two  plasma  species.  A  comparison  is  made  between  the  first 
order  correction  to  the  linear  growth  rate  obtained  in  this 
chapter  and  electron  nonlinear  Landau  growth  estimated  by 
Sloan  and  Drummond  (1970) . 

Chapter  V  deals  with  the  problem  of  electromagnetic 
waves  propagating  parallel  to  the  external  magnetic  field. 

The  theory  of  chapter  II  is  extended  to  study  four  wave  inter¬ 
actions  for  light,  whistlers  and  Alfven  modes.  A  general  for¬ 
mula  for  the  amplitude  dependent  frequency  shift  is  derived. 

As  a  special  case,  the  frequency  shift  for  the  above  waves  in 
a  cold  plasma  is  obtained. 
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CHAPTER  II 

THEORY 


2.1  Occupation  Number  Formalism  (Second  Quantization) 

In  the  quantum  mechanical  investigation  of  systems 
composed  of  large  numbers  of  particles,  it  is  beneficial 
to  use  the  occupation  number  formalism  or  'second  quanti¬ 
sation'  as  it  is  often  called.  In  this  method,  one  seeks 
a  mathematical  formalism  in  which  the  occupation  number  of 
the  states  (and  not  the  coordinates  of  the  particles)  play 
the  part  of  the  independent  variables. 


a.  Bose  Statistics 

Consider  a  system  of  N  non-interacting  particles 

•  “V  , 

(bosons)  with  wave  functions  £  (x) ,  £  (x) ,  . . . ,  which  then 

1  2 

form  a  complete  set  of  orthogonal  and  normalised  wave  func¬ 
tions.  The  total  wave  function  of  this  system  is  a  sym¬ 
metrised  sum  of  products  of  the  functions  (x) 


2 


N  !  N  ! 
1  2 

N! 


(x 


(x  ) 


PnUn) 


(1) 


Here  N.  are  the  number  of  particles  in  the  ith  state, 

l 

p  ,p  , ...p  are  the  ordinal  numbers  of  the  states  in  which 
1  2  N 

the  individual  particles  are,  and  the  sum  is  taken  over  all 


' or-  k  ^ 


- 


7 


permutations  of  those  suffixes,  which  are  different. 


constant  factor  is  chosen  so  that  the  function 

N  ,N 

1  2 


is  normalized. 


The 


Let  f  be  the  operator  of  some  physical  quantity 

Ut 

acting  only  on  a  function  of  x^.  Consider  the  operator 

f  =  y  f  (2) 

L  a 
a 


which  is  symmetrical  with  respect  to  all  the  particles  and 

whose  matrix  elements  with  respect  to  the  wave  functions  (1) 

will  be  determined.  These  matrix  elements  will  be  different 

from  zero  for  transitions  which  leave  the  numbers  N  ,N  , . . . , 

1  2 

unchanged  corresponding  to  the  diagonal  elements,  or,  since 
each  of  the  operators  f  acts  only  on  one  function  in  the 
product  (1) ,  its  matrix  element  can  be  different  from  zero 
for  transition  whereby  the  state  of  a  single  particle  is 
changed,  which  implies  that  the  number  of  particles  in  one 
state  is  diminished  by  unity,  while  the  number  in  another 
state  is  correspondingly  increased.  Simple  calculations 
yield  the  non-diagonal  elements 


N.  ,N.  -1 
l  k 


CF) 


N.-1,N, 
l  k 


= 


(3) 


where  f.,  is  the  matrix  element 
lk 


(x)  f  £k  (x) d3x 


(4) 
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The  diagonal  matrix  elements  of  F  are  the  mean  values  of  the 

quantity  F  in  the  states  XT  and  are  given  by 

N  N  .  .  . 

1  2 


f  =  y  f . .n. 

v  ill 


(5) 


The  operators  F  can  be  pictured  as  acting  on  the  occupation 
numbers  if  one  introduces  operators  a^ ,  which  decrease  by 
one  the  number  of  particles  in  the  ith  state  and  possess  the 
matrix  elements 

N.-l 

i 

(a.)  =  /N7  (6) 

N. 

l 


The  Hermitian  conjugate  operators  a*  obviously  have  the  ma¬ 
trix  elements 


(7) 


i.e.  they  increase  the  number  of  particles  by  one.  It  is 
easily  shown  that  the  operator  F  can  be  written  as 


F  =i„fikaiak 

i/k 


(8) 


The  matrix  elements  of  this  operator  are  the  same  as  those 
of  (3)  and  in  fact  this  is  the  expression  for  F  in  the  sec- 
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ond  quantisation  form.  Thus,  one  had  been  able  to  express  an 
ordinary  operator  (of  the  form  (2))  acting  on  functions  of 
the  coordinates,  in  the  form  of  another  operator  acting  on 
functions  of  new  variables,  the  occupation  number  INK.  This 
result  is  easily  generalized  to  operators  pertaining  to  more 
than  one  particle  at  once. 

From  (6)  and  (7) ,  one  can  prove  the  following  com¬ 
mutation  relations 


[a. 


+ 

lk 


=  a 


+  +  _  r 

a,  —  a,  a  .  —  o  . , 
l  k  k  l  lk 


=  0 


(9) 


The  operators  a^  and  a^  are  known  as  annihilation  and  crea¬ 
tion  operators. 

Finally,  it  remains  to  express,  in  terms  of  the  oper¬ 
ators  a^  and  at,  the  Hamiltonian  H  of  the  physical  system  of 
N  identical  interacting  particles  that  is  actually  being 
considered.  In  the  nonrelativistic  approximation  H  can  be 
represented  in  general  form  as  follows 


E  =  I  H  +  l  U(x  ,x  )  +  l  U(x  ,x  ,x  )  +  ...  (10) 

a  01  a>e  01  P  p  ' 

Here  H  is  the  part  of  the  Hamiltonian  which  depends  on  the 
coordinates  of  the  ath  particle  only  H  -  -  +  u^xa^ ' 

where  U(x  )  is  the  potential  energy  of  a  single  particle  in 

a 


' 


10 


an  external  field.  The  remaining  terms  in  (10)  correspond  to 
the  mutual  interaction  energy  of  the  particles.  In  the  sec¬ 
ond  quantisation  form,  equation  (10)  becomes 


H 


l 

i  ,k 


H 


ik 


+ 

a .  a 

l 


k 


+ 


ik£m 


[U] 


ik  +  + 
n  a  .  a.  a  n  a 
Jim  i  k  Z  m 


+  .  .  . 


(11) 


This  gives  the  required  expression  for  the  Hamiltonian  in  the 
form  of  an  operator  acting  on  functions  of  the  occupation 
numbers.  For  a  system  of  non-interacting  particles,  equation 
(11)  becomes 

H  =  £  Hxkatak  (12) 

l  ,k 


If  the  function  (x)  is  taken  to  be  the  eigenfunction  of  the 
Hamiltonian  of  an  individual  particle,  the  matrix  H^k  is 
diagonal,  and  its  diagonal  elements  are  the  eigenvalues  E^  of 
the  energy  of  the  particle.  Thus  equation  (12)  yields 


E  =  7  E.N.  (13) 

V  11 
1 

This  is  the  expression  for  the  energy  levels  of  the  system. 

The  formalism  which  has  been  developed  can  be  put  in 
a  somewhat  more  compact  form  by  introducing  the  operators 


(14) 


With  the  properties  of  a^  and  a^. 


it  is  clear  that  the  oper- 
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ator  V (x)  decreases  the  total  number  of  particles  in  the 
system  by  one,  while  XF+  (x)  increases  it  by  unity.  One  can 
easily  prove  that  these  operators  satisfy  the  commutation 
rules 

I'Hx)  .^(x')]  =  6  (x-x 1 ) 


[¥(x),nx')]  =  [V+(x)  ,f+(x')]  =  0 


(15) 


The  expression  (8)  for  the  operator  (2)  can  be  written, 
using  the  operators (14 ) ,  in  the  form 


F 


¥+  (x)  f  T  (x)  d  3x 


(16) 


which  is  the  same  as  (8).  Thus,  the  Hamiltonian  H  can  be 
written,  using  (16)  and  similar  expressions,  as  follows 


H 


d  3x¥+  (x)  H  (x)  ¥  (x)  + 


xd  ax ' (x) ¥+  (x '  )U(x,x'  )  ¥  (x)  ¥  (x '  ) 


T  •  •  •  • 


(17) 


To  clarify  equation  (17) ,  suppose  each  particle  of  the  given 
system  is  described  (at  a  given  instant)  by  the  same  wave 


function  £ (x)  which  is  normalised  so  that  — 


|5(x) |  d3x  =  1. 
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Then,  it  is  immediately  evident  that,  if  one  replaces  the 


operators  ¥(x)  in  (17)  by  £ (x) ,  this  expression  leads  to  the 


mean  energy  of  the  system  in  the  state  considered.  This 
gives  the  following  rule  for  deriving  the  Hamiltonian  in  the 
second  quantisation  formalism.  The  expression  for  the  mean 
energy  is  written  in  terms  of  the  wave  function  of  an  indi¬ 
vidual  particle,  and  this  function  is  then  replaced  by  the 
operators  ¥ (x) . 

b.  Fermi  Statistics 

The  basic  theory  of  the  method  of  second  quantisation 
remains  unchanged  for  systems  of  identical  particles  obeying 
Fermi  statistics.  However  the  actual  formulae  for  the  matrix 
elements  of  quantities  and  for  the  operators  a^  are  different, 
because  unlike  the  case  of  Bose  statistics  the  wave  functions 
are  antisymmetric. 

The  matrix  element  of  an  operator  F  of  the  type  (2) 
are  in  the  present  case 


(18) 


for  the  diagonal  elements,  and 


£(i+l,k-l) 


(19) 


for  the  off-diagonal  elements.  The  factor  (-1)  appears  to  a 


. 

L  f  m  4  H  Olfi  I 
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power  equal  to  the  sum  of  the  occupation  numbers  of  all 
states  between  the  ith  and  the  kth  (i<k) ,  £  (i  +l,k-l)  = 
k-1 

I  n  and  0. ,1.  refer  to  n.  =  0,  n.  =  1  (n.  being  the  num- 
r=i+l  r  1  1  1  1  1 

ber  of  fermions  in  the  ith  state) . 

In  order  to  represent  the  operator  F  in  the  form  (8) , 

the  operators  a^  which  will  be  replaced  by  for  fermions 

must  be  defined  as  matrices  whose  elements  are 

(c.)  1  =  (C.  )  1  =  (-1)  (20) 

h  °i 

These  operators  satisfy  the  following  anticommutation  rules: 


C.C 

l 


+ 

k 


+  etc. 

k  i 


{Ci'Ck}  =  {Ci'Ck}  =  0 


(21) 


From  the  definition  (20)  ,  it  is  obvious  that  the  result  of 
the  action  of  the  operators  depends  not  only  on  the  num¬ 
ber  n.  itself  (as  in  the  case  of  Bose  statistics),  but  also 

l 

on  the  occupation  numbers  of  all  the  preceding  states. 

With  these  properties  of  the  operators  and  C+ ,  the 
formulae  (8)  and  (11)  remain  valid.  The  formulae  (16)  and 
(17)  which  express  the  operators  of  physical  quantities  in 
terms  of  the  operators  y(x)  and  ¥+(x)  defined  by  (14)  also 
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hold  good.  The  commutation  rules  (15)  ,  however  changes  into 
the  anticommutation  relations 


{y (x) ,¥+ (x' ) }  =  6 (x-x' ) 


{ ¥  (x)  ,  Y  (x  ’ )  }  =  {Y+  (x)  ,Y+  (x’  )  }  =  0 


(22) 


2.2  Motion  of  a  Charged  Particle  in  a  Uniform  Magnetic  Field 
The  non-relativistic  Hamiltonian  for  a  particle  of 
species  j  in  a  uniform  magnetic  field  is  given  by 


H  .  =  _  .  , 

j  2m.  [  ]  c 


e  . 


P_;+77^-  A(x)  |  +  e^4)(x) 


(23) 


where  is  the  canonical  momentum  of  the  jth  particle  and 
A  (x)  and  <}>  (x)  are  the  vector  and  scalar  potentials  of  the 
field  which  can  be  split  into  time  independent  zeroth  order 
term  and  a  time  dependent  perturbation  term  which  in  this 


case  will  be  due  to  collective  excitation  of  waves. 


A  (x)  =  A 

0 

->• 

+  A 

l 

(x) 

(24) 

4>(x)  =  (p 

+  <P 

(x) 

(25) 

o  1 


Here  A  is  the  vector  potential  due  to  a  uniform  external 
o 

magnetic  field  (0  =  0);  A  (x)  and  (p  (x)  are  the  vector  and 

o  1  i 


1 
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scalar  potential  due  to  collective  fields  of  transverse  and 
longitudinal  waves  respectively. 

Substituting  (24)  and  (25)  into  (23)  yields 


where 


H.  =  Hf0)  +  H?1* 
3  3  3 


(2) 


H  : 
3 


(°)  _ 


2m 


(?j  + 


+  h: 
3 


4  a  ) 2 

c  o 


(26) 


(26. a) 


H 


(O 


e  . 


+  A  (x)  •  (P.+  —Jam  +e .  4>  (x) 
i  3  c  o  j  3  i 


(26. b) 


H  : 
3 


(2) 


e .  2 

— ^ t  (x)  , 

2m. c 2  i 


(26. c) 


and  H j ° ^  is  the  unperturbed  Hamiltonian.  The  first  and  sec¬ 
ond  terms  in  equation  (26. b)  are  responsible  for  the  linear 
interaction  of  particles  of  species  j  with  the  electromag¬ 
netic  waves  and  the  third  term  is  responsible  for  the  linear 

•  .  ( 2 ) 
interaction  of  particles  with  longitudinal  fields.  , 

equation  (26. c) ,  will  be  responsible  for  the  nonlinear  wave- 

particle  interaction  (electromagnetic  waves) . 

The  vector  and  scalar  potentials  A  (x)  and  <J>  (x)  can 

i  i 

be  expanded  in  Fourier  series  in  a  large  box  of  volume  V, 
assuming  the  usual  periodic  boundary  conditions  (Harris,  1969) 


. 


1 
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A(x)  =  l 

q 


2irhc2 


VwT(q)FT(q) 


1/2 


e  \A  e 

ql  q 


i  [q*x-ajT  (q)  t] 


+  -itq-x-u)  (q)t] 
+  A  e 

q 


(27) 


(p  (x)  =  l 
i  . 


47Tftto  (q) 


Vq2FL(q) 


1/2 


a  e 

.  q 


i [q*x-wL (q) t] 


-l  [q •  x-o)  (q)  t] 
+  a  e  L 

q 


(28) 


The  reason  for  the  factors  in  square  brackets  is  discussed 
later.  The  functions  F^ (q)  and  FL(q)  are  given  by 


FT(q)  = 


1  3  r  2  t 


w=:wT(q)  ' 


(29) 


fl(5>  = 


3w 


(u>e  (<o,q)) 


cu=0JL  (q )  ' 


(30) 


w  (q)  ,  co  (q)  are  the  longitudinal  and  transverse  frequencies 
•Li  1 

• 

respectively;  e  (q,co)  ,  e_  (q,w)  are  called  the  longitudinal 
and  transverse  dielectric  functions  and  a  ,  A  are  Fourier 

q  q 

coefficients;  a*  and  A+  are  their  complex  conjugates.  Fi¬ 
nally,  e^  is  a  polarization  vector. 

The  transition  from  classical  to  quantum  mechanics 


(quantisation  of  electromagnetic  fields)  is  made  by  reinter- 
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preting  and  as  annihilation  operators  for  longitudinal 
and  transverse  fields  (bosons)  of  momentum  ftq  and  energy  fico  ? 

q 

a*  and  A+  are  the  corresponding  creation  operators.  These 
operators  satisfy  the  commutation  relations  (9) . 

The  energy  in  the  electromagnetic  field  (transverse 
waves)  is  given 


(31) 


This  is  not  the  total  energy  associated  with  the  wave;  in  so 
far  as  the  particles  of  the  medium  move  in  response  to  the 
wave,  their  energy  must  be  properly  included  in  the  total 
energy.  Landau  and  Lifschitz  (1960)  have  shown  that  in  such 
a  dielectric  medium  the  total  energy  is 


(32) 


This  shows  that  the  energy  density  which  a  wave  would  have 


in  a  vacuum  must  be  corrected  by  the  factor 


oo2eT  (o>) ) 


(jd=0)r 


=  FT(q)  when  it  moves  in  a  medium  of  a  dielectric  constant 
e^vO)).  Substituting  equation  (27)  into  (32)  and  performing 
the  integration,  one  obtains 


UT  =  I  fiooT(q)A+Aq 

q 


(33) 


' 
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Similarly 


ul  =  h  [ 


r  , 

2 

J  d8x 

Et 

L 

^(“EL(“n 


00=00, 


=  ^  jd!x|v<H 


2  3  ,  t  \  \ 

5S(“ei(u)) 


(34) 


00  =  00, 


is  the  total  energy  in  the  electric  field  of  longitudinal 
waves.  Substituting  (28)  into  (34)  and  performing  the  inte¬ 
gration  yields 


(35) 


It  is  obvious  that  the  terms  in  the  square  brackets  of  equa¬ 
tions  (27)  and  (28)  provide  the  correct  form  for  equations 
(33)  and  (35). 

2.3  Particle-Quasiparticle  Interaction 

According  to  the  theory  developed  in  sec.  1,  the  total 
interaction  Hamiltonian  between  particles  of  species  j  and 
longitudinal  quasiparticles  can  be  written  in  the  second 
quantisation  formalism  as 


Hint  =  jd3x'i'+  (x)  ej4)1  (36) 

Here  ¥ (x)  and  ¥+ (x)  are  the  field  operators  which  satisfy 
the  anticommutation  relations  (22)  and  <J>  (x)  is  the  electro¬ 
static  potential  given  by  (28).  The  field  operators  ¥ ' s  can 

be  expanded  in  the  annihilation  and  creation  operators  C_^ 

k 


' 
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and  C+  which  satisfy  (21) 
k 


T(x)  =  l  C  £  (x) 

■+  k  k 
k 


V+(x)  =  l  C*g*(£) 

t  k  k 
k 


(37) 

(38) 


where  E,  (x)  is  the  solution  of 
k 


the  Schrodinger  equation 


(H  -  E  )  £  (x)  =  0 
0  t  it 


(39) 


and  H 

o 

plasma 


is  the  unperturbed  Hamiltonian.  In  unmagnetized 

(A  =0)  the  unperturbed  Hamiltonian  is  given  by 
o 


H 

o 


h2 

V2  and  (39)  yields 

j 


r  /"N  1  ik*X  /An\ 

(x)  =  —  e  (40) 

£  /V 

Substituting  (28),  (37),  (38)  and  (40)  into  (36)  and  carrying 

out  the  integration,  one  obtains 


Hint  =  I  M<q)C-  +  M*(q)C  _>c%± 

*  +  0  k+q  k  q  0  k+q  k  4 
k ,  q 


where 


M  (q)  =  M  (q)  =  e. 
o  o  J 


4 tt hwL  (q)  */2 


Vq2FT (q) 


(41) 


(42) 


The  terms  in  the  interaction  Hamiltonian  (41)  can  be 
represented  by  Feynman  diagrams  as  follows:  the  first  term 
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which  contains  C*  describes  the  process  in  which  a 

k+q  k  q 

particle  in  the  state  it  is  destroyed,  one  is  created  in  the 

state  k+q  and  a  quasiparticle  of  momentum  hq  is  destroyed. 

The  term  containing  C_^  ^C^a^  describes  the  inverse  process. 

k+q  k  q 

These  processes  are  represented  by  the  diagrams  in  Fig.  1  and 
yield  the  so  called  'linear  Landau  damping'. 


Fig.  1  Linear  Landau  damping 

The  interaction  Hamiltonian  for  a  magnetized  plasma 
undergoing  electromagnetic  wave-wave  and  wave-particle  inter¬ 
actions  and  its  representation  in  Feynman  diagrams  will  be 
developed  later  in  chapter  V. 

2.4  Scattering  of  Particles 

Consider  a  test  particle  of  charge  e^  and  mass  m^ 
interacting  with  a  plasma,  the  interaction  Hamiltonian  is 


ter 
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given  by 


e  .  e 


H 


int 


(X)  =  II  =  V  (x) 


(43) 


S  1  x-x  . 

SI 


—r 

where  V(x)  is  the  electrostatic  potential  energy  at  the  posi- 
tion  x  of  the  test  particle,  and  xg^  is  the  position  of  the 
ith  particle  of  species  s  in  the  plasma.  Only  Coulomb  inter¬ 
actions  are  taken  into  account.  The  potential  V(x)  can  be 
expanded  in  a  Fourier  series  in  a  box  of  volume  V  as  follows 


V  (x)  = 


l  V(q)e 

q 


-* 

iq«x 


(44) 


where  V(q)  is  given  by 


V  (q )  = 


d3x  -iq*x 
-y-  e  ^  V (x) 


=  I  I 


47re.e  -iq*x 
_ ^  51 


si  Vq 


(45) 


Thus,  substituting  (44)  and  (45)  into  (43)  yields 


4Tre  .e  iq*  (x-x,, .  ) 

I  I  I  - ^  e 


-*  S  1 

q 


Vq 


(46) 


The  total  interaction  Hamiltonian  in  the  second 
quantization  formalism  is  given  by 


H.  . 
int 


d3xd3x,'l'+(x)'t,+  (x’)Hint(x)'f  (x)H'(x')  (47) 
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Substituting  the  field  operators  (37)  ,  (38)  and  the  inter¬ 
action  Hamiltonian  (46)  into  (47)  and  performing  the  inte¬ 
grations,  one  obtains 


H 


int 


x*  4"  4* 

I  V  C  C  C  - 

,q  k  k'  k’-q  k+q 


where  the  coupling  constant  is  given  by 


V 

q 


47re  .e 


s 


(48) 


(49) 


The  interaction  Hamiltonian  (48)  ,  can  be  represented  by 
Feynman  graphs  as  follows:  a  particle  in  the  state  k+q  is 
scattered  into  a  state  5  by  a  particle  initially  in  the 
state  k'-q  which  is  scattered  into  a  state  k'.  Momentum  is 
conserved  in  the  process.  This  is  shown  in  Fig.  2(a). 


Fig.  2(a)  Coulomb  interaction 


Fig.  2(a)  is  simply  a  scattering  of  a  particle  of  species  j 
with  a  particle  of  species  s.  In  a  many  body  problem,  this 


69V19MOO 
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process  represents  only  the  first  term  in  an  infinite  series 
which  represent  the  total  scattering  that  takes  place.  If 
the  sum  of  the  infinite  series  is  represented  by  V  ^  where 
the  first  term  is  V  ,  one  obtains  in  the  random  phase  approx¬ 
imation 

Veff  =  +  y-o-  H — 


where 


4Tre  e  . 
_ s_j _ 

Vq2  e (q,w) 


£  (q,w)  =  1  -  - 


C 


(50) 


(51) 


is  obviously  the  dielectric  function  of  the  plasma  medium. 
The  effective  field  (50)  is  represented  diagrammatically  in 
Fig.  2(b). 


Fig.  2(b)  Effective  Coulomb  interaction 


■ 
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2.5  The  Scattering  Matrix 

The  theory  of  the  scattering  matrix  (S-matrix)  can  be 
developed  by  introducing  the  interaction  picture  (I.P.)  in 
which  the  state  vector  Y ( t)  is  given  by 

ifi  =  HintW(t)  (52) 

Thus,  in  the  I.P.  the  state  vectors  have  the  time-dependence 
of  Schrodinger  picture  with  the  interaction  Hamiltonian  H^nt 
instead  of  the  total  Hamiltonian.  It  is  useful  to  obtain  a 
solution  of  this  equation  particularly  suited  to  describing 
scattering  processes,  with  the  initial  state  of  the  system 
at  t  =  -  00 

y(-oo)  =  =  |  i>  (53) 

This  state  vector  will  completely  specify  the  particles 
present  initially  (i.e.  long  before  scattering  occurs  when 
all  particles  are  still  far  apart).  Equation  (52)  then  tells 
how  the  state  vector  (53)  changes  with  time.  In  particular, 
it  gives  the  final  state  ¥(+0°)  =  =  |  f>  at  time  t  =  °°, 

long  after  scattering  is  over  and  all  particles  are  far  apart 
again.  The  S-matrix  is  defined  as  the  operator  which  trans¬ 
forms  V  (-00)  into  ¥  (+°°)  : 


Y  (+oo )  =  SV  (-“) 


(54) 


■ 
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so  that  finding  the  S-matrix  is  equivalent  to  solving  equa¬ 
tion  (52)  which  can  be  transformed  into  an  integral  equa¬ 
tion  as 

w  (t )  =  ¥(-»)  +  j  dt]Hint(tj)'nti)  (55) 

—  oo 


Solving  this  equation  by  iteration,  one  obtains  the  S-matrix 


00 


n 


s  =  I 

n=0 


(-i) 

i-fcn 

n  !ti 


+  oo  +oo 

1  / 


dt  ...dt  P{H(t  ) H  (t  )...H  (t  )} 
i  n  i  i  i  2  in 


—  00 


(56) 


Here  the  Dyson  chronological  product  of  n  factors 

means  that  the  factors  are  not  to  be  taken  in  the  order  in 
which  they  are  written  in  the  P-bracket;  but,  operators 
with  later  times  stand  to  the  left  of  operators  with  earlier 
times.  The  equivalence  of  equation  (56)  holds  for  each  term 
n  separately  and  is  easily  verified. 

The  S-matrix  expansion  (56)  will  effect  many  compli¬ 
cated  transitions.  However,  only  certain  terms  will  contrib¬ 
ute  to  a  given  transition  |  i>  |f>  as  they  must  contain 
just  the  right  absorption  operators  to  destroy  the  particles 
present  in  | i>  and  that  they  must  contain  just  the  right 
creation  operators  to  emit  the  particles  present  in  |f>. 

It  is  convenient  to  introduce  what  is  called  the  nor¬ 
mal  product  (N-product)  and  the  Wick's  chronological  product 
(T-product) .  An  operator  which  is  a  product  of  creation  and 
annihilation  operators  is  called  a  normal  product  if  all 
creation  operators  stand  to  the  left  of  all  annihilation 


,  ,  ,.,v  „i  O  IX  ,'jon  l3f  K...J  si  rt6=  ri;-rtw  -  ' 


' 
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operators.  For  example 


N 


_  1  1  2  2 


+  + 

—  _  \|/  \y  \J/  \y 

12  12 


'F+H'+4'  y 
12  2  1 


(57) 


The  minus  sign  resulting  from  the  interchange  T  «--►  y+ .  Such 

1  2 

an  operator  first  absorbs  a  certain  number  of  particles  and 
then  emits  some  particles.  It  does  not  cause  emission  and 
reabsorption  of  intermediate  particles.  Thus,  one  writes 
each  term  in  the  S-matrix  expansion  (56)  as  a  sum  of  normal 
products.  Each  of  these  will  effect  a  particular  transition 
|  i>  -*■  |f>,  which  can  be  represented  by  Feynman  graphs. 

The  T-product  is  defined  as  follows: 

(i)  For  two  boson  fields,  the  T-product  is  equal  to  the  P- 
product  defined  by 


P{<j>  (x  )  (p  (x  )  } 
1  2 


4>  (x  )  <J>  (x  )  ,  if  t  >t 
12  12 


4>  (x  )  (x  )  ,  if  t  >t 
2  1  2  1 


(58) 


(ii)  If  A (x  )  is  one  of  the  operators  y(x  )  or  (x  )  and 

i  i  i 

B ( x  )  is  one  of  the  operators  T (x  )  or  V+ (x  )  then  the  P- 
2  2  2 

product  is  defined  as  (58)  but  the  T-product  is  defined  by 


T{A(x  ) B (x  ) } 
1  2 


A (x  ) B (x  ),  if  t  >t 
12  12 


-B (x  )A(x  ),  if  t  >t 
2  1  2  1 


(59) 


Because  of  the  anticommuting  nature  of  fermion  operators, 
the  T-product,  which  allows  for  the  change  of  sign  when 
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interchanging  two  fermion  fields, is  more  appropriate. 

Thus,  equation  (59)  can  be  written  in  general  as 


T{A(x  )  B  ( x  )  .  .  .  }  = 
1  2 


(-1) PP{A (x  ) B (x  ) . . . } 
1  2 


(60) 


where  the  exponent  P  is  the  number  of  interchanges  of  pairs 

of  fermion  fields  required  to  change  from  {A(x  )B(x  )...} 

1  2 

to  the  chronological  order. 

For  theories  in  which  the  interaction  Hamiltonian  is 
bilinear  in  the  fermion  field,  one  can  simply  replace  the 
P-product  in  the  S-matrix  expansion  (56)  by  T-products, 
since  only  even  numbers  of  interchanges  of  pairs  of  fermion 


fields  are  involved.  Thus,  equation  (56)  becomes 


S 


1 


n=0 


n 


(-i) 

n!hn 


.dt  T  H_(t  ) H  ( t  ) . . . H  (t  ) 

n  I  I  i  I  2  In 


(61) 


From  equation  (57) ,  it  is  obvious  that  the  N-product 
of  any  set  of  operators  has  zero  expectation  value  in  the 
vacuum  state.  For  this  reason,  one  can  introduce  the  con¬ 
traction  or  chronological  pairing  of  two  operators  as  fol¬ 
lows  : 


A  (x  )  B  ( x  ) 
1  2 


t|a(xi)b(x2)| 


N 


A  (x  )  B  (x  ) 
1  2 


(62) 


Then  the  expectation  value  of  the  T-product  in  the  vacuum 


state  will  be 


- 

' 
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A~(x  )&(x  )  =  <0|THA(x  )  B  (x  )H0>  (63) 

12  {  l  2  j 

This  is  just  a  function  that  picks  out  cases  of  non¬ 
vanishing  commutators  or  anticommutators  between  the  oper¬ 
ators  in  the  product. 

All  one  needs  now  is  a  general  procedure  for  reducing 
any  more  complicated  T-products  to  N-products  and  pairings. 
This  comes  from  Wick's  theorem  which  states  that  "any  chro¬ 
nological  product  is  equal  to  the  sum  of  all  possible  N- 
products  that  can  be  formed  with  all  possible  pairings." 

This  can  be  seen  from  the  following  example 

T{ABCD>  =  N[ABCD|  +  N  [abCD  |  +  n|  ABCbj 


+  N 

[abcd] 

+  N 

[abcd] 

+  n| 

[abcdJ 

+  n| 

[abcd] 

+  N 

La15bJ 

+  n| 

[a§cd] 

+  n[SbSdJ  (64) 

The  generalization  of  this  equation  has  been  proven  by 
Wick  (1950).  Thus,  one  can  write  the  mixed  T-products  oc- 
curing  in  the  S-matrix  (61)  as  a  sum  of  N-products.  Each 
of  these  products  corresponds  to  a  definite  process  char¬ 
acterized  by  the  operators  which  are  not  contracted.  In 
the  next  section  the  contractions ,  as  interpreted  in  terms 
of  virtual  particles  in  intermediate  states,  will  be  discussed. 


. 
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2 . 6  Feynman  Graphs 

The  expansion  of  the  S-matrix  into  normal  products  is 
quite  complex,  even  for  terms  of  low  order  in  n.  The  inter¬ 
pretations  of  these  various  terms  are  greatly  facilitated  by 
the  use  of  Feynman  graphs,  introduced  in  the  previous  sec¬ 
tions  (figs.  1  and  2) .  A  qualitative  discussion  of  the  basic 
ideas  of  Feynman  graphs  for  the  second  order  S-matrix  is  given 
in  this  section.  For  n=2,  equation  (61)  yields 


(2)  =  (-i): 

2!fr 


dt  dt  T 
1  2 


int (ti)Hint (t2J 


=  (-i)2 

2  I  fi 2 


i 


d4x  d4x  T^H.  ^  (x  )H.  ^  (x  ) 


mt 


int 


(65) 


Here  x  =  (x,t) ,  d4x  =  dtd3x  and  Hi  (x)  is  the  Hamiltonian 
density.  Substituting  H^nt(x)  =  (x)  e  j <t>  (x)  ¥  (x)  into  (65) 

one  obtains 


2  _  2 


(2 )  _  (-j)  e 
2  Ifi1 


[fd‘‘x  d4x  t/i>+(x  )  <(>  (x  )Y(x  )¥+(x  )  <(.  (x  )T(x  ) 

J  J  1  2  V  1  1  1  2  2  2  J 


(66) 


Here  e<|)  is  the  electrostatic  potential  energy,  ¥ '  s  are  the 
field  operators  and  the  chronological  T-product  is  given  by 


t/ vl/+d)  Y  ¥+d>  '^■==n!'F  4>  y  <|>  'F 

1  -  r  -  -  -  -  -  J  1  1  1  2  2  2_ 


1112  2  2 


+  N 


+  * 

V  6  Y  Y  <J>  Y 
1112  2  2 


+  N 

y+4)  'Fl'+<|>  v 

+  N 

’  +  -P 

y  <j>  y  ^  4)  y 

+  N 

y  4>  y  i  <j>  y 

1  ‘  1  1  2  2  2. 

L  1  1  1  2  2  2] 

L  1  1  1  2  2  2J 

.  .  ^ 
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r  +1 — p^-.+  ’  i 

+  n  f(|  y  r<j>  y  +  n 

L  1  1  1  2  2  2J 


- F =r+ - » 

y  <t>  y  y  <{>  y 
,1112  2  2 


+  N 


Y  <p  Y  Y 

.  1  1  1  2  2  2. 


(67) 

where  Yi  e  Y(x±)  and  <p±  =  <p>(x±).  Equation  (66)  can  be  written 
as  follows 


(2)  =  s<2>  +  s<2>  +  s<2>  +  s<2)  +  s(2)  +  s<2) 

a  b  c  d  e  f 


(68) 


where 


(2 )  (— i) 2e2  f [  4  4  r  +  +  1 

'  =  — ■  d  x  d  x  n  y%  y  y %  y 

a  21ft2  J  J  1  2  L  1  1  1  2  2  2J 


(68. a) 


(2)  _  (-i)2e 


2^2  r 


2  !ft; 


d4x  d4x  Uy%  Y  Y+d>  Y 
1  2  v  L  1  1  1  2  2 


1 

2J 


+  N 


Y+<{>  Y~~Y+<j)  Y 
.  1  1  1  2  2  2. 


(68. b) 


(2)  _  (-i)2e2 


2 !  ft  2 


d4x  d4x  N 
1  2 
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Y 4>  Y  y  <t>  y 

.  1  1  1  2  2  2j 


(68. c) 


(2)  _  (-i)  2 e2 


2  !fi: 


d4x  d4x  •( N 
1  2 


|n  ’ 


rn+TF - 

f  (ji  f  Y  <f)  Y 
1112  2  2 


r  +  * — *+  ’ 

+  NY  <p  Y  Y  <f>  Y 
L  1  1  1  2  2  2^ 


(68. d) 


(2)  _  (-i)2e2 


2HV 


riT - f=fr - 1  n 

d4 x  d4x  N  Y  (p  Y  Y  4>  Y 
1  2  L  1  1  1  2  2  2J 


(68. e) 


2  „  2 


(2)  _  (-i)ze 


2  !fi: 


r:+ 


r~T 


=~y— 1  - 

d4x  d4x  N  |  Y"r4>  Y  Y  <p  Y 
1  2  L  1  1  1  2  2  2 


(68. f) 


and  contractions  over  equal  times  operators  have  been  omit¬ 
ted  as  well  as  contractions  which  obviously  vanished  such  as 
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r  *  r+-  *  1  m  t - 1 

V  ¥  =  H  =  H  =  *T  <f>  =  Y  <b 

12  1  2  12  12  21 


r+^ 

Y  <J>  =  0 

2  1 


(69) 


As  shown  in  Fig.  1  the  interaction  Hamiltonian  H. 

int 

is  represented  by  a  graph  consisting  of  a  vertex  with  two 

fermion  lines  and  one  boson  line.  The  second  order  terms 
( 2 )  . 

S  will  be  described  by  graphs  involving  two  vertices, 

suitably  joined  together.  Those  transitions  which  a  partic- 
(  2  ) 

ular  term  can  effect,  depend  on  the  external  fields 

(that  is,  those  which  are  not  contracted)  in  this  operator, 

for  it  is  those  external  fields  which  are  responsible  for 

absorbing  initial  and  creating  final  particles. 

(  2  ^ 

The  term  S  ‘  ,  equation  (68. a)  involves  only  external 
a 

lines  corresponding  to  two  processes  of  the  types  in  fig.  1 
operating  independently  of  each  other  and  has  no  physical 
meaning . 

/  2  \ 

The  operators  1  describes  Compton  scattering  as 
well  as  other  processes.  Compton  scattering  is  represented 
by  Feynman  graphs  as  shown  in  fig.  3(a)  and  (b) .  The  con¬ 
traction  is  described  by  the  line  from  the  vertex  t  to 

1  2  i 

t  .  For  t  <t  it  describes  the  emission  of  a  virtual  fermion 
2  12 

at  t  ,  its  propagation  to  t  and  its  absorption  there.  Thus, 
1  2 

^  4  ^  f 

y  ¥  is  called  a  fermion  propagator,  and  it  describes  a  vir- 
1  2 

tual  intermediate  state.  The  other  two  processes  describe 
the  creation  and  annihilation  of  two  bosons.  This  is  known 
as  nonlinear  Landau  damping  (growth)  and  respresented  dia- 


. 
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grammatically  by  fig.  4(a)  and  (b) .  The  factor  2!  in  the 
denominator  of  equation  (68. b)  will  be  cancelled  due  to 


(a)  (b) 

Fig.  3  Compton  scattering 

the  equality  of  the  two  terms.  In  equation  (68. c) there  are 
four  external  fermion  lines.  The  bosons  occur  only  in  the 

t  t  t 


12  1 


(a)  (b) 

Fig.  4  Nonlinear  Landau  damping  (growth) 

j - — i 

contraction  4)  4>  which  is  interpreted  as  a  boson  propagator . 

1  2 

The  Feynman  graph  for  this  process  is  shown  in  fig.  5  and  is 
known  as  particle-particle  scattering. 


. 
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Equation  (68.d)  is  represented  diagrammatically  in 
fig.  6.  This  leads  to  the  self  energy  of  the  fermion.  Sim- 


Fig .  5  Particle-particle  scattering 

ilar  to  (68.b) ,  the  two  terms  are  equal  by  changing  x 
in  the  second  term. 


x 


2 


1  2 


Fig.  6  Fermion  self-energy 

Equation  (68. e)  represented  in  fig.  7,  similarly 
describes  a  boson  self  energy,  that  is  a  modification  of 


Fig.  7  Boson  self-energy 
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boson  energy  due  to  its  interaction  with  the  fermion  field. 
This  interaction  enables  the  boson  to  create  a  virtual 
fermion  pair  which  subsequently  absorbed  again. 


Fig.  8 

Finally,  fig.  8  shows  the  graph  representing  equation 
( 6 8  - f ) -  This  has  no  external  lines  at  all  and  consequently 
does  not  cause  any  transitions. 

This  completes  the  analysis  of  the  second  order  term 

( 2  ) 

S  into  normal  products  and  the  interpretation  in  terms  of 
Feynman  graphs.  The  extension  to  a  decomposition  into  higher 
order  contributions  although,  of  course,  more  complicated, 
presents  no  essential  difficulty  and  is  described  in  the  rest 
of  this  section. 


The  Feynman  graphs  for  three  and  four  wave  interac¬ 


tions  can  be  obtained  from  the  third  and  fourth  order  terms 
of  the  S-matrix.  For  n  =  3,  equation  (61)  gives 


3 !  h 

where 


S<3)  =  -til 


(70) 


e 


(71) 


« 
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and  the  terms  which  yield  the  three  wave  interactions  are 


T  { 


l^Oe  3  <  N 


r+ . ~  + . 1  +  '  "I  r  r2  -J=Y+  V - ' 

+NY  <J>  Y  ^4)  ¥  ^  cb  Y 

.1112223331  |111222333 


(72) 


where  the  factor  6  is  due  to  the  summation  over  all 
possible  fully  contracted  products.  Equation  (72)  can  be 
written  in  the  following  form  after  taking  the  N-product  of 
the  contractions 


T  { 


r  qr~:  rn  rrn 

6e  3  V  Y  ¥  ¥  y  T 

3  1  2  2  3  3 


'+  r+ : 

\p  \j/  \j/ '  \|/  \|/ 

1  3  3  2  2  1 


N'f  cj)  <p  <p 

l  1  2 


(72. a) 


Substituting  (72. a)  into  (70)  yields 


(3)_ 


(-i) 3e? 

—  -  . -  -jrL 

ft3 


d  x  d4x 
1  2 


r  r  v ''  •'X~i  ,-t~i  rp  rr~ i  -]  f 

I  n  I  +ff  ^  ^  ¥  y  ]  Nl  4>  (p  <p 

L  1  2  2  3  3  1  1  3  3  2  2  lj  {  1  2  3 

(73) 


In  equation  (73)  there  are  three  internal  fermion  lines 
(three  propagator)  and  three  external  bosons.  The  Feynman 
graphs  for  these  three  wave  interactions  are  shown  in  fig.  9. 


Fig.  9  Three  wave  interactions 


Similarly,  for  n  -  4,  equation  (61)  yields 
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where  the  terms  in  the  T-product  which  are  responsible  for 
the  four  wave  interactions  are 


T  { 


T  ’ '+  1  ’  +  *  ‘  +  ~‘  ‘  '  •  4-  f  1  4-  •  •  4-  :  ‘  + 

I  u/  U/  VI/  i"  —  * . '  • 


}  =*>  (-1 )  2  4  e  *!  ¥  m  V  y  y  y  y  T  +  y  y  Y  V  Y  ¥  ¥  y 
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1  4  4  2  2  3  3  1 
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<P  <P  <P  4>  I 

1  2  3  4  J 


(75) 


In  this  process,  there  are  four  fermion  propagators  and  four 
external  bosons.  As  shown  in  equation  (75)  there  are  six 
possible  nonequivalent  diagrams.  Each  of  these  has  4!  equiv¬ 
alent  diagrams.  Equation  (75)  is  represented  diagrammatically 
in  fig.  10.  The  choice  of  the  direction  of  the  bosons  (figs. 

9  and  10)  depends  on  the  problem  under  investigation. 

2.7  General  Rules  For  the  Calculations  of  the  Scattering 
Matrix  Element 

The  scattering  matrix  element  for  the  three  wave  in¬ 
teractions  will  be  calculated  as  an  example  from  which  the 
general  rules  for  more  complicated  diagrams  are  established. 

The  scattering  matrix  for  this  process  is  given  by  (73)  and 

rp 

the  fermion  Propagator  Y  'F  is  calculated  as  follows: 


’ 
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Fig.  10  Four  wave  interaction 
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V+V  =  < 0 | T j 1  I 0>  =  0  if  t  >t 
12  l  i  2j  12 


■<0  |  V  ¥  0> 

2  1 


if  t  <t 
1  2 


1  r 

-  "  v  *  e 


—  iE  T  -*  \ 

t  -ik* (x  -x  ) 


1  2 


(76) 


i  (k*x-E  t) 


"I  -A-  \  A  U  U  / 

where  the  field  operator  Y(x)  =  —  [  C  e  k  has  been 


used  and  T  =  t  -t  . 

2  1 


/V  *  k 
k 


At  this  point,  one  needs  an  analytical  expression  for 
a  function  with  the  following  properties  (Ziman,  1969) 


g(Ej  = 

k 


-iE  T 

\r 

-e  T>  0 


(77) 


0 


T<  0 


This  can  be  constructed  by  a  contour  integration 


+°° 


g  (e  ) 
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27ii 


dEe 


-iET 


—  oo 


E-E  +i6 
-> 

k 


(78) 


where  the  integral  is  along  the  real  axis  in  the  plane  of 
the  complex  variable  E,  and  the  small  positive  quantity  6 
is  simply  a  convergence  factor.  When  (78)  is  substituted 
into  (76) ,  it  yields  the  following  formula  for  the  fermion 
propagator 
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G(E,E+) 
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E-ErH-ifi 

k 


(79. a) 


.  .  +  +  ' 
Similarly,  one  can  write  expressions  for  ¥  V  and  ¥  ¥  . 

2  3  2  1 

For  the  specific  example  under  investigation,  one  is 


interested  only  in  the  contribution  of  (73)  for  the  specified 


initial  and  final  state  configuration.  This  is  the  annihila¬ 


tion  of  two  bosons  at  x  ,  x  and  the  creation  of  another  at 

1  2 

x  (fig.  9).  Hence,  from  the  factor  <f>  (p  (p  one  chooses  the 

3  12  3 

term  that  has  this  property.  Since  <p  and  (p  annihilate 

1  2 

bosons  and  d)+  creates  another,  one  retains  only  <t>  cj)  where 

3  12  3 

<J>  =4)  +  <{>+  ,  etc.  Therefore,  the  N-product  in  (73)  yields 
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(q  ) M  (q  )a+  a+  a-* 

2  0  3  q  q  <3 

3  1  2 


-i 


•x  -or*  t  )  -  fq 
3  q  3J  K  l 

3 


e 


(80) 


where  M  (q)  is  given  by  (42) . 
o 

r 

Substituting  the  similar  expressions  for  Y  ¥  (n,m  = 
1,2,3  and  n  ^  m) ,  from  (79)  and  (80)  into  (73),  one  obtains 


*  3 
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(81) 


The  6-functions  imply  the  conservation  of  momentum  and 
energy  parameters  at  each  vertex.  Alternatively,  (81)  gives 

,  "V  -V  t  t 

after  performing  the  k  ,k  ,E  and  E  integrations, 

2  3  2  3 
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(82) 


where  again  the  6— functions  ensure  the  conservation  of  mo 
mentum  and  energy  of  the  initial  and  final  states. 


' 


The  scattering  matrix  element  is  given  by 
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(85) 


have  been  used. 

From  the  preceding  calculations,  one  can  easily  es¬ 
tablish  general  rules  to  obtain  the  scattering  matrix  ele¬ 
ment  (83)  for  the  three  wave  interactions. 

1.  Construct  all  possible  nonequivalent  diagrams  (fig.  11) 

.  •  •  •  ->  ->  ,  “► 
with  initial  q  ,q  and  final  q  plasmons. 

12  3 

2.  Conserve  momentum  and  energy  parameters  at  each  vertex. 

The  matrix  element  will  be  the  multiplication  of 

a.  (-—)  for  each  vertex  and  (-1)  for  each  fermion  loop, 
n 

b.  a  factor  M  (q)  for  each  vertex, 

o 


3. 
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c.  a  fermion  propagator  iG(E,E  )  =  i/(E-E  +i6)  for 

k  k 

each  intermediate  state  (internal  line) , 

d.  6  (E^-E^)  due  to  the  energy  conservation  of  initial 
and  final  state, 

and  the  sum  over  the  internal  variable  £  and  the  integra¬ 
tion  over  the  energy  parameter  E. 


Fig.  11  Three  wave  interaction  example:  an¬ 
nihilation  of  two  plasmons  and  creation  of 
another. 

These  rules  are  general  for  any  diagrams  which  can  include 
any  number  of  external  waves. 

One  can  generalize  these  rules  to  include  virtual 
bosons  by  writing  a  boson  propagator  similar  to  that  of 
(3.c)  (iD(a),w  )  =  i/w-o)  +i 6)  if  the  internal  lines  consist 
of  bosons. 

The  relation  between  the  scattering  matrix  element 
and  the  transition  matrix  element  Mfi  is  given  by  (Davydov,  1969) 


■ 
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(86) 


from  (83)  and  (86) ,  one  obtains  after  performing  the  E- 
integration 
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(87) 


To  obtain  a  common  factor  f (k) ,  one  makes  the  following 

v  y  y 

transformation  k+q  -*k  in  the  2nd  term,  k-q  -*k  in  the  third 

2  1 

term,  etc.  then  after  taking  the  classical  limit 
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l  f(£)-v[f(v) 
k  1 


hk‘ 

2m 


and 


hk 


m .  v 
3 


one  obtains  upon  expansion  of  the  denominators  up  to  the 
order  of  ti2  in  equation  (87) 
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(88) 


This  is  the  same  result  obtained  by  Harris  (1969)  .  These 
calculations  will  be  used  later  in  chapter  IV  to  find  the 
nonlinear  dielectric  function  for  the  ion  acoustic  waves. 
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CHAPTER  III 

LINEAR  AND  NONLINEAR  COLLISIONAL 
DAMPING  OF  PLASMA  WAVES 


3.1  Introduction 

The  problem  of  electron  oscillations  in  a  plasma 
neglecting  interparticle  collisions  has  been  studied  amongst 
others  by  Landau  (1946)  ,  Bohm  and  Gross  (1949)  ,  Van  Kampen 
(1955)  and  Bernstein,  Green  and  Kruskal  (1957).  The  quali¬ 
tative  importance  of  collisional  damping  was  pointed  out  by 
Bohm  and  Gross,  where  they  noted  that  since  Landau  damping 
was  negligible  for  q  (wave  vector) 0,  collisional  effects 
must  dominate  in  this  limit.  Bhatnagar,  Gross  and  Krook 
(1954)  have  shown  that  in  the  presence  of  collisions,  plasma 
oscillations  decay  with  a  rate  proportional  to  the  colli¬ 
sional  frequency,  however,  a  complete  calculation  of  this 
effect  was  not  carried  out  until  Ichikawa  (1960)  examined 
the  effect  of  binary  and  tertiary  corrections  in  high  tem¬ 
perature  plasma  on  the  basis  of  the  so-called  BBGKY  equation 
for  a  system  of  charged  particles.  Similar  calculations 
based  on  the  same  equation  have  been  carried  by  Willis 
(1962).  A  formula  for  the  dielectric  function  including  the 
effects  of  the  binary  Coulomb  collisions  has  been  derived  by 
many  authors  during  the  sixties.  Among  these,  Comisar 
(1962)  investigated  the  attenuation  of  longitudinal  plasma 
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oscillations  by  transforming  the  Landau  (Fokker-Planck) 
equation  in  time,  space  and  velocity  to  obtain  an  integral 
equation  for  the  transformation  of  the  distribution  func¬ 
tion.  Ogaswara  (1963)  used  the  Boltzman  collision  integral 
and  expanded  the  distribution  function  in  a  complete  set  of 
velocity  dependent  polynomials.  Buti  and  Jain  (1965)  have 
used  Comisar's  method  to  study  the  damping  of  longitudinal 
as  well  as  transverse  plasma  oscillations.  Thourson  and 
Lewis  (1965)  calculated  the  high  frequency  conductivity  of 
a  fully  ionized  plasma  by  using  the  BBGKY  hierarchy.  They 
assumed  that  the  terms  representing  Coulomb  potential  can  be 
replaced  by  a  cutoff  Coulomb  potential  with  a  range  equal  to 
the  Debye  length. 

More  recently  this  problem  has  been  treated  by 
Shkarofsky  (1968),  Matsuda  (1969),  and  McBride  (1969) 
through  the  use  of  the  Fokker-Planck  equation  or  the  BBGKY 
hierarchy.  Shkarofsky  studied  the  problem  of  longitudinal 
as  well  as  transverse  plasma  waves  by  using  the  Fokker- 
Planck  equation  and  also  used  the  Boltzman  equation  to  study 
the  electron-neutral  collision.  Matsuda  obtained  the  dis¬ 
persion  relation  including  collisional  effects  by  applying 
the  BBGKY  equation  and  the  pair  correlation  function  given 
by  Rostoker's  supperposition  principle,  (Rostoker's  test 
particle  method) .  McBride  solved  the  linearized  kinetic 
equation  for  the  perturbed  electron  distribution  function 
for  the  case  of  longitudinal  plasma  oscillations  along  an 
external  magnetic  field,  where  the  collision  term  is  speci- 
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fied  to  be  the  Landau  (Fokker-Planck) collision  integral. 

The  problem  of  collisional  damping  has  also  been 
studied  using  quantum  mechanical  methods.  DuBois ,  Gilini- 
sky  and  Kivelson  (1963)  calculated  the  collisional  damping 
for  longitudinal  and  transverse  waves  from  the  conductivity 
which  is  derived  by  using  the  diagrammatic  techniques  of 
field  theory.  Similar  work  has  been  done  for  transverse 
waves  by  Perel  and  Eliashberg  (1962) ,  Ron  and  Tzoar  (1963) 
and  DuBois  and  Gilinisky  (1964).  For  a  basic  understanding 
of  this  technique  one  may  refer  to  Chappel  and  Britten 
(1966)  ,  Rauscher  (1968)  ,  Harris  (1969)  as  well  as  chapter 
II  of  this  thesis. 

In  the  present  chapter,  the  quantum  mechanical  ap¬ 
proach  is  adopted  and  the  collisional  damping  as  arising  out 
of  the  scattering  of  two  particles  through  effective  Coulomb 
field  and  emission  (absorption)  of  one  plasmon  is  calcula¬ 
ted  directly  without  reference  to  the  dielectric  function  or 
the  conductivity.  In  addition,  the  scattering  of  a  particle 
pair  with  the  emission  (absorption)  of  two  plasmons  will  be 
examined.  This  process  may  be  called  a  nonlinear  colli¬ 
sional  process  similar  to  that  of  nonlinear  Landau  damping 
which  involves  two  plasmons  (fig.  4).  In  section  3.2,  the 
theory  of  chapter  II  is  extended  to  calculate  the  scattering 
matrix  elements  for  the  one  and  two  wave  processes.  Sec¬ 
tion  3.3  deals  with  the  linear  collisional  damping  of  lon¬ 
gitudinal  oscillations  and  section  3.4  with  the  two  waves 
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collisional  process. 

3.2  Scattering  Matrix  for  Linear  and  Nonlinear  Collision 

The  Feynman  graphs  for  processes  involving  scattering 
of  two  particles  through  an  effective  Coulomb  field  and 
emission  of  one  plasmon  are  obtained  from  the  following  S- 
matrix 


(1) 


(x  )  H 


where  H  and  H  are  the  interaction  Hamiltonian  densi- 
w-p  p-p 


ties  for  wave-particle  and  particle-particle  scattering 
respectively.  These  are  given  by 


H 


—  I  j  ^  »  \u  \  jn.  /  x  \ 

w-p  i  3  i  i 


=  ¥  H  (x  )  e  .  c}>  (x  )  y  (x  ) 


(2. a) 


H 


p-p  J  2  3  2332  3 


(  2  .  b) 


Substituting  (2. a)  and  (2.b)  into  (1)  yields 


s (2 )  =  (~i) 2 
2  !h2 


u  a  a  j.  ]  i  c:.y  1  1  1 

1  2  3  V  1  H  1  1  2  3 


(3) 


The  T-product  can  be  written  as  follows 


■  - 
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(4) 


where  the  factor  2  is  due  to  the  equal  contraction  terms 
Substituting  (4)  into  (3) ,  one  obtains 


(O  _  (-i) 


-  d4x  d4x  d4x 

•ft  2  J  1  2  3 


vF+y  V+V+V(2 ,3)  ¥  ¥ 
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+  Similar  terms  2*-*3  ]  e  .  (J) 

^  0 


(5) 


where  the  first  term  can  be  represented  by  Feynman  graphs 
according  to  the  theory  of  chapter  II  as  follows:  at  the 
point  x  a  fermion  is  annihilated  and  another  is  created  due 

l 
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to  the  emission  (absorption)  of  a  boson.  The  created  fermion 
propagates  to  the  point  x  ,  and  is  then  scattered  into  its 

3 

final  state  due  to  the  Coulomb  collision  field  V(2,3).  It 

x 


Fig.  12(a)  Linear  collisional  process 

can  also  be  seen  that  one  fermion  is  annihilated  and  another 

is  created  at  the  point  x  .  These  processes  are  shown  in 

2 

f ig .  12(a). 

x 


Fig.  12(b)  Linear  collisional  process 

Similarly,  the  second  term  is  represented  by  fig.  12(b) 

r+-‘ 

where  the  fermion  propagator  in  this  case  is  given  by  1 3  1  • 

Since  figs. 12  (a)  and  (b)  are  nonequivalent,  their  sum  will 
yield  the  correct  S-matrix  for  the  linear  process.  These 
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diagrams  will  be  constructed  in  the  momentum  space  and  the 
calculation  of  the  S-matrix  will  be  carried  on  in  sec.  3.3. 

A  diagram  similar  to  that  of  Landau  damping  (fig.  1) 
is  obtained  from  the  third  term  of  (5).  This  is  shown  in 
fig.  13.  However  this  diagram  is  nonconsistent,  because  the 


Fig.  13  Indirect  linear  Landau  damping 


Coulomb  potential  becomes  singular  at  the  plasmon  frequency. 
The  similar  terms  in  (5)  yield  diagrams  similar  to  that  of 
figs.  12  and  13  where  the  plasmon  is  emitted  (absorbed)  by 
a  different  species  of  fermions. 

A  process  of  pair  particle  scattering  via  the  Coulomb 
field,  with  the  emission  of  two  plasmons,  can  be  obtained 
from  the  third  order  term  of  the  S-matrix  which  contains  two 
wave-particle  interaction  Hamiltonians  and  only  one  particle- 
particle  interaction  Hamiltonian.  This  is  given  by 


(6) 


As  has  been  discussed  in  the  previous  chapter,  the 
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T-product  contains  connected  as  well  as  non  significant  dis¬ 
connected  diagrams.  Each  of  the  connected  diagrams  has  3!  =6 
equivalence.  In  addition  there  will  be  several  nonequivalent 
diagrams  that  must  be  added.  Keeping  only  the  connected  dia¬ 
grams,  (6)  yields 
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(7) 

The  first  six  terms  describe  the  nonlinear  collisional 

processes.  These  can  be  represented  diagrammatically  by  six 

nonequivalent  diagrams.  Only  the  process  describing  the  first 

x 


Fig.  14  Nonlinear  collisional  process 
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term  is  shown  in  fig.  14.  The  six  diagrams  will  be  con¬ 
structed  later  in  section  3.4  The  indirect  nonlinear  scat¬ 
tering  is  described  by  the  last  two  terms  and  is  shown  dia- 
grammatically  in  fig.  15.  The  similar  terms  will  yield 


Fig.  15  Indirect  nonlinear  scattering 


diagrams  similar  to  that  of  figs.  14  and  15  where  the  plas- 
mons  are  emmitted  by  a  different  species  of  particles. 

3.3  Linear  Collisional  Damping  of  Plasma  Oscillations 

The  scattering  matrix  for  the  linear  process  is  easily 
calculated  by  transforming  diagrams  11(a)  and  (b)  to  the 
momentum  space.  In  doing  so,  one  has  to  bear  in  mind  the 
conservation  of  momentum  and  energy  of  the  initial  and  final 
states  as  well  as  the  conservation  of  momentum  and  energy 
parameters  at  each  vertex.  This  is  shown  in  fig.  16. 
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(b) 
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k 


Fig.  16  Linear  collisional  process:  The 
upper  line  refers  to  a  heavy  ion  and  the 
lower  one  to  an  electron. 


There  are  four  diagrams  like  the  one  shown  in  fig. 
16(a).  Of  these,  half  will  be  those  in  which  the  plasmon  will 
be  emitted  by  the  heavy  ion.  Such  diagrams  give  a  negligible 
contribution  (of  the  order  m/irr  )  ;  thus  only  contributions  due 
to  diagrams  which  emit  plasmons  at  the  electron  line  will  be 
considered.  Following  the  S-matrix  rules  of  chapter  II,  one 
can  immediately  write  the  S-matrix  element  for  the  processes 
shown  in  fig.  16  as  follows 
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where  the  energy  parameters  E  and  E'  are  known  and  there  is 
no  internal  variable  over  which  to  sum  or  integrate.  The 
factor  (2tt)  is  because  the  number  of  vertices  exceed  the 
number  of  internal  lines  by  unity.  From  the  conservation  of 
energy  parameters  at  each  vertex,  one  obtains 
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Putting  hk  =  mv  and  expanding  (9)  in  power  of  ii ,  one  obtains 
by  virtue  of  the  S^,  M  ^  relationship  ((86)  of  Chapter  II) 
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where  higher  powers  of  h  have  been  neglected,  since  these 
will  not  contribute  in  the  classical  limit. 

An  equation  for  the  rate  of  change  of  N->,  the  number 
of  plasmons  of  momentum  hq  can  be  written  schematically  as 
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->  -> 


This  means  that  one  adds  all  the  processes  in  which  a  plasmon 

of  momentum  tiq  is  emitted  and  subtracts  all  those  in  which 

one  is  absorbed.  The  difference  gives  the  net  increase  in 

N-*.  Taking  f  (k)  ,  the  number  of  particles  of  momentum  ftk; 

4 

the  schematic  equation  can  be  converted  into  a  mathematical 
equation  by  replacing  each  diagram  by  the  transition  proba¬ 
bility  per  unit  time  for  the  process:  Thus,  (11)  yields 
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where  M .  is  given  by  (10)  and  E+  is  the  kinetic  energy  of 
f  l  k 

the  particle  expressed  in  units  of  h.  Expanding  the  f's  and 
E's  to  the  lowest  order  of  h,  (12)  becomes; 
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Taking  N  >>  1  and  substituting 
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into  (13),  one  obtains  in  the  classical  limit 
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where  y  ,  is  the  rate  at  which  the  wave  amplitude  decays, 

IjU 

q.q'  =  qq'cos9,V  and  v'  are  the  electron  and  ion  velocities 
respectively.  The  volume  of  the  box  in  which  the  system  is 


. 


58 


quantized  tends  to  infinity  so  that  the  sums  go  over  into 
integrals,  hence  (14. d)  and  (14. e).  Finally,  f  (v)  and  f (v ' ) 
are  the  classical  particle  distribution  functions.  If  elec¬ 
trons  and  ions  are  in  thermal  equilibrium  with  Maxwellian 
distribution,  (15)  becomes 
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Performing  first  the  v ' -integration ,  then  using  the 
asymptotic  expansion  method  (Fried  and  Conte,  1961)  to  per¬ 
form  the  v-integration ,  one  obtains 
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which  yields  after  performing  the  q'  and  0  integrations 
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Substituting  00^  =  (o^e 
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into  (18)  yields 
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This  result  agrees  with  previous  calculations  (DuBois,  Gil- 
inisky  and  Kivelson  (1963),  Shkarofsky  (1968),  Matsuda 
(1969,,  and  McBride  (1969)). 

3.4  Nonlinear  Collisional  Process 

In  the  previous  section,  the  Coulomb  scattering  of  a 
particle  pair  via  emission  of  one  plasmon  is  investigated. 

It  has  been  shown  that  for  waves  of  long  wave-length  (q+0) , 
the  collisional  damping  dominates  over  Landau  damping.  How¬ 
ever,  it  appears  that  the  nonlinear  collisional  process  is 
also  more  important  than  the  nonlinear  Landau  damping  in  the 
long  wave  length  limit. 

In  studying  the  nonlinear  collisional  process,  it  can 
be  shown  that  similar  to  the  linear  collisional  process 
described  in  section  3.3,  diagrams  involving  emission  of 
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plasmons  by  the  heavy  ions  will  make  a  negligible  contribu¬ 
tion,  and  also  there  are  six  nonequivalent  diagrams  which 
will  give  a  significant  contribution.  These  diagrams  are 
described  by  the  first  six  terms  of  (7)  which  in  the  momen¬ 
tum  space  yields  the  diagrams  shown  in  fig.  17. 

Following  the  rules  of  the  scattering  matrix  explained 
in  chapter  II,  and  the  previous  section,  one  can  easily  write 
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Fig.  17  Nonlinear  collisional  process:  The 


upper  line  refers  to  a  heavy  ion  and  the 
lower  one  to  an  electron. 
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Substituting  E  = 


fik 


'k  2m-"'  =  mv  ;'-nto  (20)  ,  expanding 

the  denominators  up  to  the  order  of  ft2 ,  and  using  relation 

(86)  of  chapter  II,  one  obtains 
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In  obtaining  (21),  terms  of  order  ft0,  and  ft  vanished 
identically  and  that  of  order  ft2  yielded  the  correct  clas¬ 
sical  limit  of  the  matrix  element.  Inserting  the  square  of 


(21)  into  the  rate  equation  for  N  yields 
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which  can  be  rewritten  after  converting  summations  into 
integrations  and  taking  the  classical  limit 
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where  I(q^)  -  tiaHj  N-*  is  the  plasmon  energy  in  the  mode  q.  . 

*i  q^ 

Considering  the  special  case  of  longitudinal  waves  where 

a)-)-  -  w->  -(0  ,  and  assuming  that  the  electrons  and  ions 

q  q  pe  ^ 

1  2 

obey  a  Maxwellian  distributions  as  before,  (23)  yields  after 
performing  the  velocity  integrations 
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where  Y^c  t^ie  nonlinear  collisional  damping  (growth)  , 
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and  0  is  the  angle  between  q+q^  and  q'  as  shown  in  fig.  18. 

Substituting  (25)  into  (24)  and  performing  the  q'- 
integration  one  obtains 
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The  angles  a  and  3  are  shown  in  fig.  18  and  <p ,  <J>  are  the 

2 

•  -  , 
azimuthal  angles  of  q  and  q  respectively.  In  order  to  per- 

2 

form  the  q  -integration,  the  plasmon  energy  I (q  )  will  be 

2  2 

considered  independent  of  q  ,  and  is  replaced  by  I(t),  the 

2 

average  value  of  I  (q  ) .  This  approximation  is  good  enough 

2 

for  plasmons  which  are  not  far  removed  from  equilibrium. 

#  , 

Furthermore,  the  integration  over  q  will  be  cut  off  at  q  , 

2  L) 
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the  Debye  wave  vector.  Thus,  for  an  average  energy  spectrum, 
(26)  gives 
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where  terms  of  order  q2  have  been  neglected. 

The  nonlinear  contribution  in  (28)  is  of  a  higher  order 

In  equilibrium  the  plas- 
mons  obey  Bose  distribution  giving  T  =  xT  and  rendering 


in  the  plasma  parameter  (n  A3)  1 

o  D 


Y 


NLC 


=  C.  When  the  plasmons  are  slightly  below  or  above  the 


equilibrium  level  of  energy,  the  above  equation  shows  that 
there  is  growth  or  decay  respectively.  In  either  case,  the 
nonlinear  collisional  process  tries  to  restore  equilibrium. 
This  is  quite  visible  from  the  solution  of  the  differential 
equation  (28)  which  yields 
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where  I  is  the  equilibrium  value  of  I (t)  at  t  -  0. 
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(29) 


68 


CHAPTER  IV 


A  STABILITY  MECHANISM 


FOR  ION  ACOUSTIC  WAVES 


4.1  Introduction 

The  problem  of  linearly  unstable  ion  acoustic  modes 
has  been  studied  by  Kadomtsev  (1965)  and  more  recently  by 
Sloan  and  Drummond  (1970).  The  ion  acoustic  waves  become 
unstable  when  the  component  of  electron  beam  velocity  in  the 
direction  of  wave  propagation  exceeds  the  phase  velocity 


(Jackson,  1960)  ,  where  T^  is  the  electron  tempera- 


c 


s  /  m . 
v  1 


ture  and  itk  is  the  ion  mass.  In  this  case,  direct  electron- 
ion  Coulomb  interactions  are  negligible  as  has  been  discussed 
by  Kadomtsev  and  Pogutse  (1967) .  Kadomtsev  argued  that  non¬ 
linear  Landau  damping  of  ion  acoustic  waves  by  the  ions  re¬ 
sembles  the  damping  of  Langmuir  waves  by  the  electrons,  and 
therefore,  he  suggested  ion  nonlinear  damping  as  the  mechan¬ 
ism  for  terminating  the  instability  of  the  ion  acoustic 
modes.  Thus,  the  electron  nonlinear  Landau  growth  was  thought 
to  be  of  little  significance  for  ion  acoustic  instability. 

This  result  was  later  used  by  Sagdeev  (1967)  to  calculate  the 
well-known  Sagdeev  formula  for  anomalous  resistivity.  Sloan 
and  Drummond  (1970),  questioned  the  validity  of  Kadomtsev's 
latter  assumption  and  therefore,  the  anomalous  resistivity 
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results  of  Sagdeev.  They  calculated  the  electron  nonlinear 
Landau  growth  and  compared  it  with  the  ion  nonlinear  Landau 
damping  and  concluded  that  the  dominant  nonlinear  mechanism 
is  electron  nonlinear  Landau  growth  which  would  be  further 
destabilizing.  It  therefore,  became  necessary  to  consider 
other  processes  such  as  mode  coupling  (Drummond  and  Pines, 

1962;  Vedenov  et  al.,  1962;  Aamodt  and  Drummond,  1964).  In 
mode  coupling,  however,  transfer  of  energy  from  one  mode  to 
another  may  arrest  the  growth  of  one  unstable  mode  at  the  cost 
of  exciting  others  and  is  not  necessarily  a  stabilization 
mechanism. 

Since  the  ion-acoustic  mode  plays  an  important  role  in 
the  heating  of  charged  particles,  for  example,  in  the  magne¬ 
tosphere  causing  auroral  precipitation  (Buneman,  1959;  Swift, 
1965;  Kennel  and  Petschek,  1969) ,  as  well  as  in  laboratory 
plasmas  (Jones  and  Alexeff,  1965;  Alexeff  et  al.,  1968; 

Okamoto,  Tamagawa,  1971) ,  this  phenomenon  will  be  further  in¬ 
vestigated  in  this  chapter.  As  a  stabilization  mechanism  for 
the  ion-acoustic  modes,  the  amplitude  dependent  frequency 
shift  (Fukai,  et  al.  1969,  1970)  which  is  of  the  same  order  of 
magnitude  as  the  nonlinear  Landau  damping  (growth)  will  be  con¬ 
sidered  and  its  consequences  discussed.  It  will  also  be  demon¬ 
strated  that  the  nonlinear  dielectric  properties  of  the  medium 
modifies  the  wave-particle  coupling  constant  and  that  it  has 
stabilizing  influence  on  the  wave  (Selim  and  Krishan  to  be  pub— 


, 
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lished  in  the  Phys.  of  Fluids,  1971-1972).  In  section  2, 
the  frequency  and  growth  rate  of  the  waves  for  the  linear 
approximation  is  derived  and  section  3  contains  the  calcula¬ 
tions  of  the  nonlinear  dielectric  function  and  the  amplitude 
dependent  frequency  shift.  In  section  4  a  comparison  between 
the  result  obtained  in  this  chapter  and  that  estimated  by 
Sloan  and  Drummond  (1970)  is  made. 

4.2  Linear  Dielectric  Function  and  Growth  Rate 

In  the  Feynman  diagram  representation  as  shown  in 
chapter  II,  the  bare  Coulomb  potential  as  well  as  the  ef¬ 
fective  potential  (which  takes  into  account  the  screening  of 
the  Coulomb  field)  are  given  by 


where  e  is  the  electronic  charge,  V  is  the  volume  of  the 
system,  q  is  the  wave  vector  for  the  qth  Fourier  component 
and  £  (q , oo )  is  the  linear  dielectric  function  of  the  plasma. 
In  finding  c(q,oo)  in  the  random  phase  approximation,  one 
adds  an  infinite  series  such  as 


n  ;  r-  4 
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-O 


by  using  the  Dyson  equation, 
obvious  that  the  linear 


47Te2 

Vq2 


1 - (3) 


From  equations  (2)  and  (3)  it  is 


dielectric  function  can  be  written  as 


e  (q,aj)  =  1  -  l  - * - O 

j 


(4) 


where  j  is  the  summation  over  the  plasma  species  with  ref¬ 
erence  to  electron  and  ion  loops.  The  'bubble'  in  the  series 
which  is  the  second  term  in  (4)  will  be  calculated  by  using 
the  S-matrix  theory  which  has  been  developed  in  chapter  II. 
Accordingly,  use  will  be  made  of  the  following  rules: 

(1)  in  the  expression  of  the  second  order  matrix  element  a 


factor 


47Te' 

Vq2 


appears  in  the  product  for  each  Coulomb  line, 


(2)  a  factor  iG(E,E^)  =  E_E^+f5  appears  for  each  internal 


particle  line,  where  E+  and  E  are  energy  and  energy 
parameters  of  the  particle  in  units  of  h. 
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(3)  sum  over  the  internal  variable  k  and  integrate  over  the 
internal  energy  parameters  E, 

(4)  multiply  by  (-1)  for  each  fermion  loop  and  introduce 

6  (E^-Ef)  which  appears  in  the  theory  and  serves  to  con¬ 
serve  energy  between  initial  and  final  states, 

(5)  use  the  relationship  ((86)  of  chapter  II)  between  the 
S  and  M  matrices  to  find  the  matrix  elements  of  the 
second,  third,  fourth  and  other  terms  in  (3). 

Applying  these  rules  to  the  calculation  of  e(q,oa),  one  ob¬ 
tains 
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Performing  the  E-integration  yields 
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A  considerable  simplification  can  be  made  by  the  linear 
transformation  k+q-^k  in  the  second  term  since,  in  any  event, 

— y  _ 

the  summation  is  over  all  values  of  k.  Thus,  there  is  a 
common  factor  f ^ (k)  which  is  the  number  of  particles  in  the 
state  k.  Therefore,  (6)  yields 
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by  substituting  E+ 


nk-l  and  hit  =  m.v  after  changing  it- 
2m .  J 
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summation  into  v  integration,  (7)  thus  simplifies  into 

(v) 

d3v  — t  -  ,  (8) 

(w-q • v) 2 


e  (q,a))  =  1  -  l 


4Tre' 

m  . 


The  distribution  functions  for  electrons  streaming 
with  velocity  u  and  ions  with  zero  mean  velocity  are 
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where  vg/  v^  are  the  electron  and  ion  thermal  velocities 
respectively.  Substituting  (9)  and  (10)  into  (8)  and  per¬ 
forming  the  velocity  integrations  by  asymptotic  expansion 
(Fried  and  Conte,  1961),  the  linear  dielectric  function  for 

.  t 

a  wave  of  frequency  go  and  wave  vector  q  is  given  by 
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where  A  =  (T  /4iTe2n  )!/2,  the  Debye  wave-length  and  go  .  , 
De  e  0  P1 

the  ion  plasma  frequency.  The  dielectric  function  given  by 
(11)  is  identical  to  the  results  obtained  by  standard  tech¬ 
niques.  The  frequency  go  can  be  written  as  w  =  go-*  +  iY:N 

SI  SI 

which  when  substituted  in  the  equation 


'  ♦ 


. 
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£  (q  /  w)  =  0 


(12) 


will  determine  the  frequency  of  oscillations  and  the  growth 
rate  for  the  ion  acoustic  wave.  Solving  (12)  for  co->-  and  y 

q  q 

yields 


U)-*  =  q  / —  =  qc 
q  q  J m.  qcs 


(13) 


,  c 

IT  S 

Yq  =  J 8  ~  <q‘U  "  “q> 


(14) 


where  y  >  0  for  q*u  >  o)+  corresponds  to  growth  of  the  wave. 

q  q 

This  means  that  the  ion  acoustic  mode  is  unstable  in  a  cone 

of  q-  space  centered  about  the  direction  of  the  drift  velocity 

u  with  an  apex  angle  2a  given  by  cosa  =  c  /u. 

s 

Another  approach  through  which  one  can  implement  deri¬ 
vation  of  the  growth  rate  y  ,  is  the  use  of  the  rate  equation. 
If  N(q)  is  the  number  of  phonons  in  state  q  and  f  (k)  is  the 
number  of  electrons  in  state  k,  then,  the  rate  equation  for 
the  sound  waves  can  be  written  as 


dN  (q) 
dt 


2  TT 

h2 


I 

k 


M  | 2N (q) q 
o 


a£e(g) 

3k 


5  (E^-  ->■— (^“^"Er^-) 
k+q  q  k 


=  2,  © 


d  3 v  I  M  |2 
o 


N  (q) 


->  3f  (v) 

a.  e 

m 


9v 


6  (q  •  v-co-0 

q 


(15) 


where  the  direct  matrix  element  M  (derived  in  chapter  II) 

o 

-K,  1  /  2 

is  of  the  order  (^J  .  This  yields  the  correct  classical 

Substituting  the  electron  distribution  function  (9) 


limit . 


. 
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into  (15)  and  performing  the  velocity  integration,  one  obtains 


1 

N  (q) 


dN  (q) 
dt 


rn  vi 
0 


(q*u-uHj) 

qv 


(16) 


which  is  the  same  as  equation  (14)  if  one  substitutes 


4TTe2hm-> 

1/2 

M  = 

0 

q 

II 

(u>£  (q  ,  to )  ) 

_  Vq2F+ 

^3 

(16. a) 


In  calculating  F^j,  the  linear  dielectric  function  and 
the  linear  frequency  of  the  waves  have  often  been  used.  How¬ 
ever,  for  the  case  of  unstable  oscillations,  one  should  use 
the  nonlinear  frequency  of  the  waves  and  the  nonlinear  di¬ 
electric  function  in  calculating  F-*.  These  modifications 

^3 

reduce  the  coupling  constant  and  consequently  the  growth  rate. 


4.3  Nonlinear  Dielectric  Function  and  Frequency  Shift 

In  order  to  calculate  the  nonlinear  dielectric  func¬ 
tion  e  (oj,q),  one  has  to  know  the  nonlinear  processes  (J. 

N  J_j 

Coste,  1969) .  The  basic  processes  responsible  for  the  non- 
linearities  in  the  medium  are  those  where  the  bare  Coulomb 
potentials  represented  by  dotted  lines  in  (3) ,  is  replaced 
by  the  effective  potentials  which  are  represented  by  beaded 
lines  and  the  simple  loops  in  the  same  equation  are  replaced 
by  the  square  and  triangular  loops  shown  in  figs.  19  and  20. 
Then,  following  the  method  of  the  previous  section,  one  ob- 
tains  a  formula  for  eNL(oj,q). 


' 
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Fig.  19  Nonlinear  effects  in  the  dielectric 


function  (square  loops) 


tiJ1 
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q' 


(d) 


Fig.  20  Nonlinear  effects  in  the  dielectric 
function  (triangle  loops) 
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In  the  nonlinear  limit,  the  effective  potential  VNL 
is  given  by  the  sum  of  an  infinite  series  as  shown  diagram- 
matically  in  the  following  equation 


Summing  the  infinite  series  (17) ,  one  obtains 
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47Te2 

Vg  2  e 
4  NL 


(20) 


where 


eNL(^'“) 


the  nonlinear  dielectric  function  is  given  by 


£ 


NL 


(21) 


In  the  above  equation  £^  is  the  linear  dielectric  function 
given  by  (11) , 


q' ,  j 


i 

q'  j 


a  -  Sum  of  the  diagrams  in  fig.  19/V 
i  ef  f 

(22) 


a  =  Sum  of  the  diagrams  in  fig.  20/V  , 

2  eff 


(23) 


and  the  numerators  of  (22)  and  (23)  are  the  transition  matrix 
elements  of  the  interactions  shown  in  figs.  19  and  20.  These 
matrix  elements  can  be  calculated  with  the  help  of  the  rules 

(  4  ) 

developed  in  chapter  II.  Let  S  _  ^  be  the  S-matnx  element 

corresponding  to  the  diagrams  in  fig.  19,  then  one  can  write 
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00 
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^  7T  e 2 

2 

y 

4TTe  2 

q 

fi4 

Vq2eL  (q  ,  a) ) 

l 

j  rq' 

Vq  ,2 

F+, 

q  . 

N(q')  (-1) 


x 


dE 


(E-E++i6)  (E-ar>, 

K  CT 


+  Similar  term  with  ur* 

q 


-E+_^,  +  i6)  2  (E-co-rn^, 


6  (Ei-Ef ) 


~Er>-  -*■  -*  ,  +  i  6  ) 
k-q-q 

(24) 


where  (-i/h) 4 ,  due  to  the  four  vertices  in  the  diagram; 

( 47Te2/Vq 2 e  )  ,  is  the  effective  Coulomb  field; 

1  /2 

(4ne2ha)-^ , /Vq '  2F->,  )  ,  is  the  coupling  constant  due  to  wave- 

particle  interaction;  (-1) ,  is  because  of  the  single  fermion 
loop,  the  four  energy  denominators  are  due  to  the  four  vir¬ 
tual  particle  lines  and  N(q')/  is  due  to  integration  over  the 
energy  parameter  of  the  virtual  plasmons  because,  there  are 
N(q')  plasmons  in  the  state  q' .  Thus,  the  summation  over  the 
plasmons  states  q'  simply  becomes  the  summation  over  the  wave 
vector  q' .  Finally,  the  summation  j  is  over  the  electron  and 
ion  loops,  and  the  last  term  is  due  to  the  existence  of  plas¬ 
mons  propagating  in  the  opposite  direction,  i.e.  a  plasmon 
can  be  emitted  where  it  was  absorbed  before  and  vice  versa 
without  violating  the  conservation  of  momentum  and  energy 
parameters . 

Similarly,  one  can  write  the  scattering  matrix  elements 


for  diagrams  b  and  c. 
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+  similar  terms  with  on-*, 
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6  (Ei-Ef ) 


(26) 


Performing  the  E-integration  in  (24)  one  obtains 


( )  _  4Tre 
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f j (k-q ' ) 
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(27) 


In  equation  (27),  transforming  k-q'-q  -*  k  in  the  second  term 
and  k-q'  -*  k  in  the  third  and  fourth  terms,  permits  rewriting 
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Similarly,  performing  the  E-integration  and  making  the  neces¬ 
sary  transformation  in  k,  equations  (25)  and  (26)  yield 
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hk: 


In  the  classical  limit,  one  substitutes  E.  =  _ 

k  2m . 

3 

and  changes  the  k-summation  into  integration  to  obtain 
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The  S-matrix  for  diagrams  d,  e  and  f  can  be  obtained 
from  a,  b  and  c  respectively  by  making  the  transform  q  ^  -q 
and  co  ^  -co 
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6  (Ei-Ef) , 


(34) 


(*)  _ 


—  2  77  i  >  4ire2 
^  h  j 


h 


4iTe: 


Vq2  e. 


2  fioaTr.NCq') 

l  -JL; - 

j,£’  ^  Fq' 


d  vf . (v) 
3 


,  ft  q  '  2  ^ 

w** ,  -q  *  v+x-3 — 

q  2m. 

k  3 


2  +  -*,ftq2'! 

“-q-v+2S7 

3  J 


“-s'H&r 


3 


-> ,  -*  ,  ftq  '  2 

L“$-‘q  •v+2E7- 

3 


+ 


-> ,  -*  ftq 
(!)-►,  -q  •  v — - 


i  2^  2 


t  q 


2m  . 
3 


“>■  "fl  ■>x2 

WJ.-W- (q*-q)  •v-^j-(q'-q) 

j 


3 


->  ->  ftq2  1  2  f 

“”q#v‘2Sr 

3; 


(q  -q)  *v+2FT(q  q) 


6 (Ei“Ef )  ' 


(35) 


CO  _ 


=  Pr) 


47Te2 

4Tre2 

2  ftco+,N(q') 

v  q _ 

/• 

h3 

Vq2£L. 

l  : 

.  q  '  2  F-> , 

j  /  q  q 

« 

d  3vf ^ (v) 


+  ,'ftq2 
w"q  *V+2S- 

3 


->  ,  ,  hq  ’  2  ^ 

i“5*“q  9V+W7 
14  3 


.  ->  -> 


WJ.+W- (q+q’ ) *v+2^-(q+q' ) 

q  j 


2^ 


88 


1 


±  -> ,  ftq 
oo-q  •  v+— — 


2  ^ 


2m  . 
3' 


i  -* ,  -*  ftq 
m^.-q' -v — ^ 


"\  ( 


1  ^  ->  -* 

WJ.-W-  (q'-q)  * v_2^— (q ' ~q ) 

j 


2m . 
3 


2 


+ 


1 


-*  ->  hq 2  1  f  -*■ ,  +  hq  '  2 
2m. J  (  q'  M  2m j 


q 


fS  -> 

-  (q'+q)  • v_2mT(q  +q) 


->  hq 

w-q  •  v-^-3— 
2m  . 
3 


2)  ( 


->  ,  -*■ .  hq  1 

*v+2^: 


2 


3  ' 


^,-03- (q'-q)  •v+2^-(q'-q) 
L  4  j 


6  (Ei-Ef ) 


(36) 


By  adding  the  S-matrix  elements  as  obtained  from  the  six  dia¬ 
grams  in  fig.  19  one  obtains  the  transition  matrix  element  M  . 
In  the  above  relations,  the  factor  (-2771/11)  which  appears  in 

the  S-matrix  will  be  cancelled.  Since  ha)->,N(q')  is  a  clas- 

'‘d 

sical  quantity  defining  the  wave  energy  I(q')  and  the  clas¬ 
sical  limit  of  the  M-matrix  has  to  be  independent  of  "ft,  one 
has  to  expand  the  denominators  in  equations  (31)  -  (35)  up 
to  the  order  of  ft3.  In  doing  so,  terms  of  order  of  ft0,  h  and 
-ft2  must  vanish  identically.  After  some  tedious  algebra,  this 
trivial  condition  has  been  verified  and  therefore,  the  matrix 

element  M  is  given  by 
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where 
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A  =  oo  -  q  •  v  ,  A  ' 


-  q ' • v , 


(38. a) 


and  0  is  the  angle  between  q  and  q' . 

Following  the  same  rules,  one  can  write  the  S-matrix 
elements  for  processes  shown  in  fig.  20  and  consequently  the 

transition  matrix  element  M  . 
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where  the  factors,  (4Tre2/Vq2 e  )  ,  is  due  to  the  effective  Cou¬ 
lomb  field;  [  4-ne  2fuoq  ' /Vq  '  2F+  ,  ]  1  / 2  ,  is  the  wave-particle 
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coupling  constant;  47Te2/V  |  q+q  '  |  2  e  (5+q  '  ,  w+oo  ,  )  ,  is  the  ef- 

q 

fective  Coulomb  field  line  between  the  two  triangles  and 
(-i)5/fi5/  due  to  the  five  vertices.  The  six  energy  denom¬ 
inators  are  due  to  the  six  virtual  fermion  lines,  and  the 
'similar  terms'  have  the  same  meaning  as  in  the  square  process. 
Also,  the  S-matrix  element  due  to  diagram  (b)  of  fig. 

20  is 


(-i)  5 

47Te2 

2  47re2fujo->  ,  N  (q  '  ) 

v  q  ^ 

4fTe 2 
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Vq2eL 
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dE 
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I 
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T\ 


6(E.-Ef) 


(40) 


Performing  the  E ' -integration 

the  necessary  transformation 

bk 2  ->*  — 

stituting  E^  =  ,  ftk  =  m^v 

into  integration,  one  obtains 


in  (39)  and  (40),  and  making 

-y  ->  -y  -y  _  . 

of  k-q  -q  k,  etc.,  and  sub- 

and  changing  the  k-summation 
by  adding 


s ( 5)+s  ^ 

a  b 


-2Tri^ 

(4ire2 )  2 

2  ftoo+,N(q') 

v _ q _ 

1 

h  J 

Vq2eL 
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l 

s 


h  -  2 


d3vfs(v>(^) 


s  (w-q-vj  (co+.-q'  -v)  (u)+u)+ ,  -  (q+q  '  )  *v) 

4  Si 
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(w-q • v) 


(03->,-q'  -v) 
Si 


+ 


,  -  (q+q  )  •  V 


+  . 


6  (Ei“Ef) 


(41) 


Similarly,  the  sum  of  diagrams  (c)  and  (d)  yield 


S(5)+S'5)  =  (^|i' 

c  d  v  fi 


2  \  2 


(4Tre  ) 


Vq  2  e. 
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L  I  q*j 
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[2uh 


l 


'hco^.N  (q1  ) 

1 
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+ 
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(w-q • v) 


(oj->  -q  '  •  v)  o)+co->  ,  -  (q+q  '  )  » v 

q  q 


+ 


6(ErEf) 


(42) 


By  adding  (41)  and  (42)  and  carrying  out  the  same  procedure 
as  for  the  previous  case,  one  obtains 


M 

2 


( 4TTe2 ) 


Vq  2  £  ( q  ,  (i) ) 


l 

q  Vq 


ftu)-* ,  N  (q  '  ) 

_ 2_ _ _ _ 

q  '  2  F-> , 
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1 


I 


,,,  -  ,1 
q+q  '  |  e  (q+q  '  ,oj+03+  ,  )  [  j 


I  2  + 


f.  (V)  V 

d  3  v  — J -  Y  /• 


m  . 
3 


l  I 
,1 


+ 


±  ,±  ± 


/ 

r  f  .  (v) 

2 

d3v  J  Y 

I- 

u  J 

m2  2 

3  J 

(43) 


where 


Y!  XX'  (X  +  X' ) 


-> 


q'  -  (q+q1  )q2  q-  (q+q '  ) q 1 2  q»q 1  (q+q 1 )  2 


->  -> 


X 


X 


+ 


X  +  X  ' 


Y 

2 


XX'  ( X-X ' ) 


q' • (q-q1 )q 

x 


+ 


->•  ->*  —y  n 

q  * q  (q-q  ) 

x-x ' 


(44) 


and  X,X'  are  defined  by  (38). 

The  matrix  element  M  can  be  obtained  directly  from 

2 

the  calculations  given  in  the  last  section  of  chapter  II  if 
one  recognizes  that  each  diagram  in  fig.  20  is  the  square  o 
three-wave  interaction  (apart  from  some  coupling  constant 
factors).  Substituting  (36)  and  (43)  into  (22)  and  (23) 
respectively  yield 


a  = 


1 


03 


n  Vq"£  3 
o  t 


l  ^  l 

q'^q 


i  (q 1 ) 


o  q' F£> 


d  3  vf j  (v) A j  (q / q '  ,v)  , 


(45) 


and 


a 
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V^£L  5^qq'2Fq' 
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q+q 
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1  1  V 

- — - \l 
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These  equations  occur  as  equations  (15)  and  (18)  in  the  work 
of  Selim  and  Krishan  (1971-1972). 

The  velocity  integrals  in  (45)  and  (46)  are  calculated 
by  the  asymptotic  expansion  method  (Fried  and  Conte,  1961) 
where  the  ion  integrals  converge  asymptotically  in  the  limit 
v . 

~  <<  1.  On  the  other  hand  the  electron  integrals  converge 

c 

asymptotically  in  the  limit  —  <<  1.  In  order  to  perform 

e 

,  t  ->■ 
the  velocity  integrals,  one  assumes  for  simplicity  that  q' 

t  •  . 

lies  within  a  cone  of  angle  0  <<  1  about  the  q  direction 

o 

which  is  nearly  parallel  to  the  direction  of  the  electron 
beam  velocity.  Furthermore,  the  calculations  will  be  re¬ 
stricted  to  a  spectrum  which  is  nearly  constant  within  an 
interval  about  qD  and  zero  outside,  where  qD  is  the  ion  De¬ 
bye  wave  number.  Also  it  will  be  assumed  that  the  energy  is 
of  the  form 

1(5' ,t)  =  I ( t ) 6 ( 0-0  )  (47) 

o 

This  choice  of  spectrum  is  more  general  than  that  of  Sagdeev 
and  Galeev  (1969) ,  where  they  have  considered  only  two 

.  •  i 

spectrum  lines  m  q  -space. 

Substituting  equations  (9) ,  (10)  and  the  spectrum 

distribution  (47)  into  (45)  and  (46)  and  integrating  over 
the  velocities  and  wave  vector  q' ,  one  obtains 


. 
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and 


a 


q  ADe£L 


22 4tt 2  I  ( t) 


27N 


D  mv 


>  3/2 


T . 

l' 


0 


a 


<3  ADe£L 


2  56tt  2  I(t) 


9N 


D 


mv' 


T. 

l 


3/2 


6 


(48) 


(49) 


4  TT 

where  =  —  nQ^  De  is  the  number  of  electrons  in  Debye 

sphere.  In  obtaining  the  result  for  a  and  a  ,  terms  of 

1  2 

order  higher  than  (q/qD)  have  been  neglected  and  to  keep 

(48)  and  (49)  accurate  up  to  order  I  (t)  only/  the  linear 

wave  frequency  is  substituted.  Since  0  <<1/  one  can 

o 

easily  see  that,  the  contribution  from  the  triangle  dia¬ 
grams  (49)  exceeds  the  square  diagrams  (48)  by  a  large 

factor  (1/0  )  and,  therefore,  the  former  will  be  retained. 

o 

Substituting  (49)  into  (21)  yields 


£NL 


256tt2 

9N  0 
D  o 


I(t) 

mv2 


(50) 


where  e  (q,w)  is  the  real  part  of  the  linear  dielectric 

function  given  by  equation  (11).  Solving  equation  (50)  for 

£  =  0  as  has  been  done  in  the  previous  section,  and  re- 

NL 

taining  only  terms  of  order  I (t) ,  one  obtains 


q 


qc  +  Ao)->- 
^  s  q 


where  the  frequency  shift  is  given  by 

4 


(51) 
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Au)-> 

q 


1287T2 
9N  6 

u  0 


I  (t) 
mv2 


(52) 


Next,  to  obtain  the  function  F-*  in  the  nonlinear  approxima- 

q 

tion,  one  has  to  differentiate  (21)  first  with  respect  to  oj  , 

and  then,  substituting  oj  =  qc  and  performing  the  q'-integra- 

s 

tion.  Equations  (21) ,  (45) ,  (46)  and 


F+  = 
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(“EMT  I,= 


NL 


(jO==U3- 


(53) 


yield  (appendix  A) 
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(54) 


where 


a  = 


TT2  I(t) 


0oWD  mv2 


T  . 


3/2 


(55) 


Substituting  (54)  into  (16. a)  yields  the  modified 


matrix  element 


M 


4ne2  ho)-> 

2  _ _ a 

Vq2F-> 


4l,e2R.  x2  qc 
2V  DeH  s 


,,  256  x  3/2  t  128 
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(56) 


then,  inserting  (56)  and  (51)  into  (16)  to  get 


=  /V8  qcs(/i) 
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128c 
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Thus,  since  I  (t)  increases  with  time,  one  notices 
that  the  factor  within  the  square  bracket  (which  is  obtained 
from  the  nonlinearities  associated  with  the  matrix  element) 
increases  and  consequently  the  growth  rate  decreases.  Also, 
the  expression  in  the  curly  bracket  (where  the  last  term  is 
due  to  the  frequency  shift)  decreases  as  I (t)  increases  until 
the  wave  energy  assymtotically  approaches  its  maximum  value, 


I 

max 

mv2 
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"oND 


128tt 


rT .  \ 
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T 

v  eJ 


3/2 


(58) 


making  the  growth  vanish  and  thereby  stabilizing  the  wave. 
This  is  so,  because  the  phase  velocity  of  an  initially  un- 


Fig.  21 
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stable  wave  has  shifted  from  the  unstable  part  of  the  electron 
distribution  to  the  damping  part.  This  is  shown  schematically 
in  fig.  21. 

As  a  result  of  the  instability  there  will  be  some 
heating  of  the  plasma  particles  over  any  finite  time.  To  dis¬ 
cuss  this  heating  mechanism,  one  writes  the  rate  equations 
for  electrons,  as  well  as  ions,  from  which  the  rise  in  tem¬ 
peratures  can  be  estimated  in  the  stationary  state.  The  rate 
equations  for  particles  can  be  written  as 


3£j(v)  _  2* 


9t 


—  I  lM  I  2 
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where  M  is  the  modified  matrix  element  (55)  and  is  ex- 
o  K 

pressed  in  unit  of  ft.  By  expanding  f  (kiq)  ,  Ej^+^,  and  the  6 


function  in  powers  of  ft,  and  substituting  Ej> 
(59)  yields  to  the  lowest  order  in  fi 


m .  v 
3 


1 


(60) 
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Substituting  the  value  of  and  ftw^N  (q)  =  I(q,t)  into  (60), 
and  changing  the  q- summation  into  integration,  one  obtains 
the  rate  equations  for  the  plasma  particles 


9f i (v) 
9 1 


i (q,t) 
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(62) 

Multiplying  both  sides  of  equations  (61)  and  (62)  by  mv2  and 
rruv2  respectively  and  integrating  the  result  over  velocity 
space,  one  obtains  the  rates  of  change  of  temperature 
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To  eliminate  the  first  term  in  (63) ,  multiply  (61)  by  v  and 
integrate  over  the  velocity  to  obtain 
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then,  substituting  (65)  into  (63)  yields 
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^  e 


I(q) 


(66) 


Solving  (64)  and  (66)  for  —  >>  1,  one  obtains 

c 

s 


T 

e 


1/5 


2/5 


(67) 


where  T  is  the  initial  electron  temperature, 
o 

Since  the  maximum  value  of  the  last  factor  is  unity, 
one  can  write 


T.  ra.  ’/5  c  4/5 

_i  4  (_i)  f_S) 

T  ^  lm  ;  ( u  ' 

e 


This  result  is  different  from  the  one  obtained  by  Sagdeev 
(1965)  where  it  was  shown  that 


e 


However  it  has  been  shown  by  Sloan  and  Drummond  (1970)  that 
(69)  is  incorrect. 

4.4  Comparison 

The  stabilization  theory  discussed  in  the  previous 
sections  can  be  strengthend  by  proving  that  the  amplitude 
cjgpgjKjgji t  term  in  the  growth  rate  obtained  in  the  previous 


. 
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section  is  larger  than  the  electron  nonlinear  Landau  growth 
obtained  by  Sloan  and  Drummond  (1970) .  From  equation  (5)  of 
the  above  reference,  one  obtains 


dN  (g) 
dt 


d  3 q" d  3 vN (q ) N (q"  ) A  6  (w* -q '  «v)  x 


3fe(?) 

8v 


(70) 


where  w '  =  m+oo"  ,  q'  =  q+q"  ,  A  =  ^ — — ,  and  B  is  given 

i  2  ow 

q 

by  equation  (6)  of  the  above  reference.  Substituting  the 
value  of  e(q',(jL)')  into  A,  and  calculating  B,  one  gets 
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then,  inserting  (71)  and  (72)  into  (70)  yields 
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The  number  of  plasmons  in  the  state  q"  is  given  by 
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De 


Substituting  (74)  into  (73)  and  performing  the  velocity  inte¬ 
gration,  (73)  becomes 


' 


101 


N  (q) 


=  2y 

dt  yNL 


/2 


7T 


.qc, 


n  ^  n 
o  De 


(— ) 

^  V  ; 


d3q" 

q'q" 2 


Kg") 

mv2 


(q ' •u-w' )  (75) 


Performing  the  q" -integration  by  taking,  as  before,  the  spec¬ 
trum  form  I(q" )  =  I  (t) 6  (6-0  ),  (72)  yields 

o 


where 


3  2  TT 
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From  (57),  the  dominant  nonlinear  term  in  the  growth 
rate  is  given  by 
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Comparing  (76)  with  (78) ,  one  obtains 
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fTil 

YC 
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lT 

^  eJ 
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(79) 


Thus,  the  nonlinear  term  due  to  the  electron  nonlinear  Lan 
dau  growth  is  negligible. 
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CHAPTER  V 

FREQUENCY  SHIFT  IN  ELECTROMAGNETIC  WAVES 

5.1  Introduction 

The  Krylov-Bogoliubov-Mitropolskii  (1947,  1961)  tech¬ 
niques  familiar  in  nonlinear  mechanics  have  been  used  by  many 
authors  to  study  the  behaviour  of  some  types  of  plasma  waves. 
In  particular,  Montgomery  and  Tidman  (1964)  considered  the 
nonlinear  propagation  of  an  electromagnetic  wave  in  a  cold 
plasma  and  found  an  expression  for  the  second  order  fre¬ 
quency  shift  for  both  travelling  and  standing  waves.  Tidman 
and  Stainer  (1965)  have  calculated  frequency  and  wave  num¬ 
ber  shifts  for  nonlinear  waves  in  a  finite  temperature 
plasma  considering  both  cyclotron  waves  and  electron  plasma 
oscillations.  Boyd  (1967)  obtained  an  amplitude  frequency 
shift  for  extraordinary  waves,  i.e.  with  a  propagation  vec¬ 
tor  perpendicular  to  the  direction  of  the  external  magnetic 
field.  Sluijter  and  Montgomery  (1965)  and  Das  (1968)  used 
the  same  techniques  to  obtain  the  amplitude  dependent  fre¬ 
quency  shift  for  both  light  waves  and  extraordinary  modes 
for  which  the  relativistic  corrections  are  included.  All 
the  above  authors  based  their  calculations  on  the  dynamical 
equations  for  a  cold  plasma  together  with  Maxwell's  equa- 

This  method  will  be  first  briefly  described.  Then, 


tions . 
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the  quantum  mechanical  method  will  be  employed  to  calculate 
the  frequency  shift  in  light,  whistlers  and  Alfven  waves. 

The  non-relativis tic  equation  for  momentum  conservation  is 
given  by 


3v  ->•  3v  .  e  1 

3"t  +  v * 7->  +  i<E  +  c  vAB)  =  0 
3x 


(1) 


and  the  corresponding  Maxwell's  equations  are 


V  *E  =  47re  (n  -N)  , 
o 


(2) 


VAE  + 


1  3B 
c  3 1 


=  0, 


(3) 


va$  =  -  || 

c  3 1 


4  7T 


eNv , 


(4) 


where  N  is  the  electron  number  density;  v,  the  electron 

velocity;  E,B,  the  electric  and  magnetic  fields  and  n  is 

o 

the  uniform  background  positive  ion  density.  To  the  above 
set  of  equations,  the  Krylov-Bogoliubov-Mitropolskii  (KBM) 
perturbation  expansion  is  applied 


N  =  n  +en  ^  °  ^  (a,  (f>)  +e2n  ^  ^  (a  ,  <f> )  +  .  .  .  .  (5) 

o 

v  =  0+ev  ^  °  ^  (a ,  4>)  +  e  2v  ^ 1  ^  (a  ,  <p )  + .  .  .  .  (6) 


E  =  0+eE°  (a,cf))+e22  ^ 1  *  (a,cf>)  +  -  .  .  . 


(7) 


104 


B  =  B^  +  eB^0^  (a  ,  (j) )  +e  2B  ^  1  ^  (a, $)  +  ....  (8) 

in  which  the  perturbation  correction  terms  are  functions 
of  the  wave  amplitude  a  and  the  phase  <f> ,  and  e  <<  1,  is  the 
expansion  parameter.  In  the  KBM  theory  one  uses 

3  d) 

=  w+e 2 A (a) +0  (e  3 )  ,  (9) 

|^r  =  e2D(a)+0(e3)  ,  (10) 

in  which  A  (the  frequency  shift  of  that  mode  under  investi¬ 
gation)  and  D  are  determined  by  the  requirement  that  there 
be  no  secular  terms  (terms  ^cj))  in  the  perturbation  expansion. 

In  (1969) ,  Tam  developed  a  perturbation  scheme  by 
which  a  small  amplitude  dispersion  effect  can  be  accounted 
for.  This  scheme,  based  on  the  method  of  Krylov,  Bogoliubov 
and  Mitropolskii  (1961),  takes  advantage  of  the  fact  that 
the  frequency  of  a  nonlinear  oscillator  depends  on  the  ampli¬ 
tude  of  oscillation.  Accordingly,  for  small  amplitude 
plasma  waves,  the  frequency,  wave  number  and  amplitude  can 
be  represented  by  a  power  series  in  the  energy  density  of 
the  waves  as  follows: 

ti o  =  Lo  +(jl)  (q)  a2+oo  (q)  a  +  .  .  .  .  (11) 

0  2  *+ 

where  the  first  term  is  the  linearized  result  and  the  ad- 
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ditional  terms  represent  the  effect  of  nonlinearity.  Tam 
applied  the  above  technique  to  the  investigation  of  whistler 
modes.  He  concluded  that,  if  whistlers  are  propagating 
parallel  to  the  external  magnetic  field,  the  frequency  shift 
is  zero,  and  that  the  linearized  solution  is  exact  for  the 
full  nonlinear  two-component  cold  plasma.  However,  it  has 
been  demonstrated  by  Selim  and  Krishan  (1971)  that  these 
observations  are  incorrect;  and  later  Sudan  (1971)  proved 
that  Tam's  result  is  restrictive. 

In  this  chapter,  a  general  theory  has  been  developed 
for  the  studying  of  whistlers,  Alfven  and  transverse  waves. 
The  various  physical  processes  responsible  for  the  amplitude 
dependent  frequency  shift  are  represented  by  Feynman  graphs 
using  the  theory  developed  in  chapter  II.  The  wave-particle 
and  two-wave-particle  interaction  Hamiltonians  are  obtained 
in  sections  2  and  3  respectively.  Section  4  contains  the 
calculation  of  the  scattering  matrix  element  and  a  general 
formula  for  the  amplitude  dependent  frequency  shift.  These 
shifts  for  light,  whistlers  and  Alfven  waves  are  given  in 
section  5. 

5.2  Wave-Particle  Interaction 

The  first  two  terms  given  by  (26. b)  of  chapter  II 
represent  the  electromagnetic  wave-particle  interaction. 

The  corresponding  interaction  Hamiltonian  can  be  written  in 
the  second  quantisation  formalism 
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H.  =  l 
int  L 


e  . 


2  cm . 
3  3 


f+(x) 


e  . 


(P  .+  -4  A  (x)  )  •  A  (x) 
J  c0  1 


+  A  (X)  •  (P  .+  J-  A  (x)  ) 
i  K  3  c  o  ; 


Y  (x)  d  3 x 


(12) 


Here,  A  (x)  is  the  vector  potential  due  to  uniform  external 
o 

magnetic  field  directed  along  the  z-axis  such  that  X  (x)  = 

o 

(-B  y , o , o ) ,  and  ¥ ' s  are  the  field  operators 


y(X>  =  ^  Cn ,  k  5n,K(x) 

n  ,k 


(13) 


^+(X)  =  ^  <(K(x) 

n ,  k 


(14) 


C+  and  cj  are  the  annihilation  and  creation  operators  of 
.K  f  n  k.  f  n 

•  •  ,  . 
particle  m  state  k  with  quantum  number  n,  and  E,  t  ls  t^ie 

n ,  K 

solution  of  the  Schrodinger  equation 


2m  . 
3 


P.+  - 
.  3  c 


e  . 

3  * 


A 


^n,k  En,k^n,k 


(15) 


and  is  given  by  (Landau  and  Lifshitz,  1958) 


i  (k  x+k  Z ) 

«n,k  ,k  (x)  =  L  Gn(y-yc)e 
X  z 


(16) 


The  Harmonic  oscillator  function  G  (y-y  )  has  the  form 

1 1 


/in  .go  . , 

G  (y-y  )  =  — - - 

n  c  [n!2n/2ir1>j  ^ 


77R 


Hn^jucj/n  <y-yc>l'e 


-LmjWcj/2n(y-yc)2 


(17) 
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where  y 
J  c 


ftk 


x 


m .  03  . 

3  cj 


I  e  .  B 
l_J_ 


m .  c 
3 


and  is  the  Hermite  poly 


nomial.  Inserting  (16)  into  (15),  one  obtains  the  correspond¬ 
ing  eigenvalues 


E 


n  ,k 


x 


ho)  .  (n+l/2)  + 
^  J 


(18) 


This  means  that  the  particle  behaves  like  an  oscillator  per¬ 
pendicular  to  the  magnetic  field  and  like  a  free  particle 
parallel  to  it. 

Substituting  (13)  and  (14)  into  (12)  yields  for  elec¬ 


trons 


H.  . 
int 


2mc 


l 


it ,  k '  ,n  ,n  ' 


Cit,n(t)Ck',n' 


.+ 

’k  ,n 


eB 


(P+ 


ye  )  •  A 
J  x  i 


-*■ 

A 


(P+ 


eB 


-> 

ye 


x 


■k 


■  n 


d  3x 


(19) 


Since  the  problem  under  investigation  concerns  transverse 

circularly  polarized  electromagnetic  waves  propagating  paral- 

let  to  the  external  magnetic  field  B  ,  the  vector  potential 

o 

A  given  by  equation  (28)  of  chapter  II  will  be 

l 

A  (x)  =  l 
1  -> 
q 

where,  in  (e  ±ie  )  the  plus  sign  refers  to  the  clockwise  po- 

x  y 

larization  and  the  minus  to  the  counterclockwise.  The 


27rfic 


VaJT(q)FT(q) 


1/2 


(e 


+ 


x 


ie  ) 

Y 


/2 


A  (t)eiqz+A+(t)e"iqz 


(20) 
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whistler  mode  have  a  clockwise  polarization,  while  Alfven 
and  light  waves  can  be  polarized  in  the  clockwise  or  counter¬ 
clockwise  direction.  Substituting  (16)  and  (20)  into  (19) 
and  performing  the  x  and  z  integrations,  one  obtains 

1/2 


H .  , 
int 


2mc 


I 


47Tftc2 


v^T(q)FT(q) 


I 


k  ,k  ,n,n'Ckx'kz+q'nCk'kZ'n'Aq 
X  z 


dy 


(“ce)Gn(yyc)(ryc)Gn’(ryc> 


±  hG  (y-y  )  G  ,  (y-y  ) 
n  2  c  dy  n 1  J  c 


+  H  •  C 


(21) 


To  perform  the  y-integration ,  one  can  make  use  of  the  fol¬ 
lowing  relations: 

1/2 


mw 


(-k^)  (y-yJG  ,  (y-y  )  = 


h 


c  n 


+.^Gn--r  (22) 


3?  Gn(y-yc>  =  V 


ftm 

I  ce 

ft 


r  ,/a  G 

-P1  G  , 

J  2  n-l 

'J  2  n+1 

(23) 


moj 

where  G  (y-y  )  =  G  [a (y-y  )] ,  a  =  - 

n  c  n  c 

(23)  into  (21)  yields 


ft 


.  Inserting  (22)  and 


1/2 


int 


eft 

2mc 


I 


q , k  , k  ,n  ,n  ' 

X  z 


8  7T  c  2  moo 


ce 


VcoT  (q)  Ft  (q) 


*“k  ,k  +q,n<“k  ,k  ,  n'AG 
X  Z  ^  X  z 


u 


dy/nT  Gn(y-yc)Gn,_1- 


dy/n'+i  Gn(y-yc)Gn«+i 


Using  the  integral 


-+H  •  C . 
(24) 


dyG„  (y-yc)Gn,  (y-y^)  =  « 


n 


n  ,n 


1  • 


(25) 


one  obtains 
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H .  , 
int 


eh 

2mc 


I  l 


q  k  ,k  ,n 
^  x  z 


87Tc2ma)  n 
_ ce  + 

tt  /  \  T71  /  n  k  ,k  +q,n-lCk  ,k  ,n  ^ 
Va3T(q)FT(q)  x'  z  x'  z' 


i - : - 

8ttc  moo  (n+1) 

Vu)T(q)FT(q)  '  kZ+q  'n+lCkx  '*2  ' 


+  H-C.  (26) 


In  (26)  the  first  term  is  due  to  the  clockwise  polarization, 
and  the  second  term  to  counterclockwise.  These  two  terms  are 
represented  diagramatically  in  fig.  22(a)  and  (b)  respectively. 


(a)  (b) 

Fig.  22  Electromagnetic  wave-particle  interaction 


5.3  Two-Wave-Particle  Interaction 

The  interaction  Hamiltonian  given  by  equation  (26.3) 
of  chapter  II  is  responsible  for  wave-wave-particle  inter- 
Xn  the  second  quantisation  formalism,  this  inter 
action  Hamiltonian  can  be  written 


*  ' 


110 


f 


H 


int 


Y+(X) 


2mc 


(A  «X  )  y  (x)  d3x 


(27) 


where  the  ¥ ' s  and  A^  (x)  are  given  by  (13),  (14)  and  (20)  re¬ 
spectively.  Substituting  these  equations  into  (27)  and  per¬ 
forming  the  integration  one  obtains 


H.  4- 
mt 


TTe2h  y 
mV  ^ 

q*q 


<q)“T 


(q')Fm(q)Fm(q') 

-i-  ± 


-l/2 


w  i,  -1  .  i  C  ,  n  A  A  |  tC  ,  C  ,  ■■  A  A  , 

n,kx,Kz+q  -q  n,kx,kz  q  q'  ",k  .k  -q-q'  n,kx,kz  q  q' 


+  H  •  C . 


(28) 


(28) 

term 

cess 


can  be  represented  by  Feynman  diagrams  as  follow:  The 

-+  -  ,+■ 

„  .  .  »k 

X  Z  -1  -  X  z 


containing  C  .  ,  ,  .  C  .  .  A  A  describes  the  pro- 

^n,k,k+q'-qn,k,kqq 


in  which  one  particle  in  the  state  n,k  ,k  is  annihilated, 

X  z 


Fig.  23  Two-wave-particle  interaction 
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and  one  is  created  in  the  state  n,k  ,k  +q'-q;  an  electromag- 

X  z 

netic  wave  with  momentum  hq'  is  destroyed  and  one  with  momen¬ 


tum  hq  is  created.  Similarly  the  second  term  describes  the 
process  of  destroying  one  particle  and  creating  another  with 
the  creation  of  two  electromagnetic  waves.  These  processes 
are  shown  in  fig.  23(a)  and  (b)  . 

5.4  Scattering  Amplitude  and  Frequency  Shift 


In  order  to  obtain  a  second  order  term  in  the  fre¬ 


quency  shift,  a  four  wave  interaction  must  be  considered. 
The  total  Hamiltonian  can  be  written  (Harris,  1969) 


H  =  H  +H ' 


(29) 


+  H  •  C . 


(30) 


H'  =  £  M  A  A  A  A  +  l 


++„+„+ 


M  A+A+A+A 


112  3  4 


2  12  3  4 


31234  *+1234 


+  7  M  A  A+A+A+  +  l 

L*  r  1  O  O  lx 


5  12  3  4 


6  12  3  4 


7  12  3  4 


8  12  3  4 


+  H  •  C  . 


(31) 


where  A. 

1 
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q 


4 


are  the  frequencies  of 


the  four  waves  in  the  absence  of  interactions.  The  perturbed 
frequencies  are  contained  implicitly  in  H'  which  includes  all 
the  possible  four-wave  interactions. 


Since  H'  contains  expressions  involving  four  plasmons, 


the  appropriate  Feynman  representation  must  include  diagrams 
with  four  external  plasmons.  These  are  given  by  three  dif¬ 
ferent  types.  First,  diagrams  similar  to  that  of  fig.  10  in 
which  there  should  be  always  two  plasmons  annihilated  and  two 
created.  The  conservation  of  momentum  must  be  satisfied  and 
the  quantum  numbers  obey  selection  rules  given  by  the  inter¬ 
action  Hamiltonian  (26)  at  each  vertex.  Second,  since  there 
are  creation  and  annihilation  of  two  plasmons  at  the  same 
vertex  as  shown  in  fig.  23  due  to  the  term  (A*A*) ,  it  is  pos¬ 
sible  to  construct  diagrams  similar  to  that  of  fig.  9  with 
two  plasmons  at  one  vertex  and  one  at  each  of  the  other  two 
vertices.  Third,  loops,  corresponding  to  the  second  order 
matrix  elements,  with  two  plasmons  at  each  vertex  are  also 
possible.  Retaining  only  terms  containing  the  annihilation 
of  two  plasmons  and  creation  of  another  in  (29) ,  one  obtains 


■ 


■ 
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+  1  Vq  Aq  Aq  Aq  +  H‘C- 

12  3  4 


(32) 


If  the  waves  are  taken  to  be  monochromatic  i.e.  q  =q  =q  =q 

12  3  4 

(Montgomery  and  Tidman,  1964)  ,  (32)  becomes 


H  =  +  MA+A  A+A  +  H*C 

2qqq  q  q  q  q 


IfiAqAq(“q  +  IT  AqAq>  +  H’C 


JfiAqAq(,J)q  +  +  H‘C- 


(33) 


where  Aw  is 

q 


the  frequency  shift, 


given  by 


|M  A+A 
n  q  q 


2M 

ft2oo 

q 


ftw  N (q) 
q  n 


2MI (q) 
h2oo 

q 


(34) 


and  M  is  the  sum  of  the  matrix  elements  for  the  three  possible 
processes  discussed  earlier  and  it  has  to  be  at  least  of  order 
of  ft2  . 

The  matrix  element  for  the  first  process  shown  in  fig. 
24  for  the  clockwise  polarization  can  be  written  as 


MI  = 


-  f-M 

^2lJ  1 


ft  ^  (-i) 


h4 


4  l2m.; 
3  3 


y  4  ,8  Trm  .00  M2 
eh'  '  i 


3 .  .9.3. 

Voo  F 

q  q 


(-1) 


n  (n+1 )  i  4 


l  IdE  - 

kx'kz'n  I  fE_Ek  ,n^  (E+^q"Ek  +q,n-l^E  Wq  Ek  -q,n+l^ 
L  z  z  z 


' 
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(a) 


(c) 


Fig.  24  The  contribution  to  the  scattering 
amplitude  due  to  f our-plasmon-particle  ver- 
tices.  This  comes  from  p*A  term. 
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n (n-1 )  i 


(E_Ek  ,,J  ^E+VEk  +q,n-li2(E+2VEk  +2q  ,  n-2^ 


+ 


(n+1)  z  i 


2-4 


(E"Ev  J2(E-«  ' -E. 


k  ,  n- 
z 


q  k  -q,n+l- 
z 


(35) 


where  the  factor  (-ft/2iri)  comes  from  the  relation  between  the 
S  and  M  matrix.  The  source  of  various  factors,  as  have  been 
explained  in  chapter  II.  The  three  different  terms  in  the 
bracket  are  due  to  the  three  different  diagrams  in  fig.  24. 
Performing  the  E-integration,  (35)  yields 


Mi  =  l  & 

J  J 


8Trm  .  oo 


J  cj 


Voo  F 
q  q 


I  h 

kx-kz'n 


n  (n+1) f  (k  ,k  ,n) 
x  z 


(E.  +oo  — E,  ,  }  2  (E.  w  -E.  ) 

k  k  ,n  q  k  +q,n-l;  ^  k  ,  q  k  -q,n+l; 


(E,  +oo  — E.  )  (e  —oo  — E  ,  ,  )  2 

k  k  ,n  q  k  +q,n-l;  ^  k  ,n  q  k  ~q,n+lJ 
z  z  z  z 


+ 


n  (n+1) f  (k  ,k  +q,n-l) 
X  z 


fEn  n -00  - E .  )  2  (E.  .  , -2w  -E.  ) 

1  k  +q , n- 1  q  k  ,nJ  1  k  +q,n-l  q  k  -q,n+l; 
z  z  z  z 


+ 


n  (n+1) f (k  ,k  -q ,n+l) 

X  z 


(Ek  -q,n+l+Wq  Ek  ,n^  ^Ek  -q,n+l+2aJq  Ek  +q,n-l^ 
z  z  z  z 


+ 


n  (n-1 )  f  (k  , k  , n) 
X  z 


(Ev  .n  +  a3rr  Ek  +a,n-l^  ^k  ,n+2a)q  Ek  +2q,n-2^ 


'kz,n  q  kz+q 
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n (n-1 ) f (k  ,k  +q ,n-l) 

X  CJ 


(E 


k^+q/n-1  ^q  '"‘k  »n^  ^Ek  +q,n-l'^q  ""k  +2q,n-2 


-03  -E, 


+0)  -E, 


-J 


n  (n-1) f  (k  ,k  +q,n-l) 

X  z 


K  +„  ,-r“„-Ev  J  (E 


kz+q,n-l  “'q  ,nJ  '•1Jk  +q,n-l+wq  Ek  +2q,n-2^ 


+ 


n(n-l)f(k  ,k  +2q/n-2) 
X  z 


(e  „  ~-2o3  -E  }  (E  o~03  ~E  ,  )  2 

^  k  +2q,n-2  q  k  ,nJ  v  k  +2q,n-2  q  k  +q,n-l; 
z  z  z  z 


2  (n+1)  f  (k  ,k  ,n) 

X  z 

(e,  -03  -E.  3 

v  k  ,n  q  k  -q,n+l; 
z  z 


2  (n+1) 2f  (k  ,k  -q,n+l) 
X  z 

(E.  _,,+<*)  -E.  )  3 

1  k  -q , n+ 1  q  k  ,n; 
z  z 

(36) 


After  making  the  necessary  transformation  in  k  and  n  (for 

example  k  +q  ->  k  and  n-1  n  in  the  third  term,  etc.)  , 
z  z 

(36)  can  be  rewritten 


Mi  =  l  ^ 

3  3 


8Trm.03  . 

3  C3 

V  03  F 

q  q 


i 


k  ,k  ,n 
x  z 


2ftf  (k  ,  k  ,  n )  x 
X  z 


n  (n+1) 


(E.  +03  -E 


kz,n  q  kz+q 


, n-1^  ^Ek  , n  Wq  Ek  -q,n+l^ 


n  (n+1) 


(Ek  ,n+U)q  Ek  +q,n-l^Ek  ,n  Wq  Ekz  q,n+l- 


+ 


(n+1) (n+2 ) 


(Ek  ,  n  Wq  Ek  -q ,  n+1^  ^Ek  ,n  2u)q  Ekz  2q,n+2 
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+ 


n(n-l) 


^Ek  ,n+a)q  Ek  +q,n-l^  (Ek  ,n+2aq  Ek  +2q,n~2 ^ 
z  z  z  z 


(n+1) 


n 


(E,  -03  -E  )  3 

k  ,n  q  k  -q,n+lj 


(E,  +oo  -E.  ,  ,  )  3 

1  kz / n  q  kz+q,n-l; 


ftk2  fik  (37) 

Substituting  E  =  oo  .  (n+1/2)  +  - — — ,  —  =  v  and  changing 

k  ,  n  cd  2m.  m.  z  ^  ^ 

z  3  3 

the  summation  into  integration,  (37)  becomes 
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n (n+1) 
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A  + 

3 
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3  2m, J 
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0  -J  nJ 

(38) 


where 


A  .  =  oo  +  oo  .  -  qv 
3  q  C3  z 


(38. a) 


Since  the  matrix  element  must  be  of  order  of  h2 ,  one  has  to 
expand  the  denominators  of  (38)  up  to  the  order  of  h3  to 


obtain 
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(39) 


Substituting  fio3  .n  -  ^  m.V^  and  taking  the  classical  limit, 

c  3  ^  3 

one  obtains 
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(40) 

Following  the  same  procedures,  one  can  write  the  matrix 
element  for  the  second  process  shown  in  fig.  25. 
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REk  ,JhE'VEk  -q.n+J. 


(41) 


Performing  the  E-integration  yields 
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Fig.  25  The  contribution  to  four-wave  scat¬ 
tering  from  two  plasmon-particle  vertices 
and  one  two-plasmon-particle  vertix.  This 

f?  — 

comes  from  A  and  p*A  terms. 
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II 
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X  z 
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^  k  -q,n+l  q  k  ,nJ 
z  z 


(42) 


Making  the  necessary  transformation  to  obtain  a  common  factor 

f (k  ,k  ,n) ,  one  obtains 
x  z 


M 


II 


2  2  8ttoj 

ft  V  TT6  f  e  >  CJ  y 

4  V  ^2m.J  V03  F  f  . 

3  3  q  q  kx,kz,n 


f  (k  ,k  , n ) 
x  z 
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n+1 


(e..  +03  ~E  i  )  2  (e.  ”0)  “E1  ,  )  2 

^  k  ,n  q  k  +q,n-l;  ^  k  ,n  q  k  -q,n+l; 

z  z  z  z 


n+1 


^Ek  ,n  Wq  Ek  -q,n+l^ 
z  z 


+ 


n 


(e,  +03  -E  )  2 

1  kz,n  q  kz+q,n-lj 


=  0 


(43) 

This  means  that  the  contribution  due  to  fig.  25(b)  exactly 
cancels  that  of  fig.  25(a). 

The  matrix  element  due  to  the  third  possible  process 

shown  in  fig.  26  is  given  by 
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(b) 


Fig.  26  The  contribution  to  four-wave  scat¬ 
tering  from  two-plasmon-particle  vertices . 

This  comes  from  the  A2  term  in  the  Hamiltonian. 


Performing  the  E-integration  in  (56),  it  is  obvious 
that  the  first  term  is  equal  to  zero.  This  means  that  the 
process  shown  in  fig.  26(a)  gives  zero  contribution.  Thus, 
only  process  26(b)  contributes  to  ^ ,  which  becomes 
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+ 


f (kx/kz+2q'n) 


(E1  ,  _  —  2 oo  -E_  ) 

v  k  +2q,n  q  k  ,n; 
2  z 


(45) 


Making  the  transformation  k  +  2q  -*  k  in  the  second  term  and 

z 


substituting  E^  ,n  =  go__.  (n+1/2)  + 


CD 


ftk2 

■x — -,  hk  =  m.v,  into  (45) 
2m .  z  j 


and  changing  summation  into  integration,  one  obtains 
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tt  e 2 

2 

‘  h2  J 

1 

h->0  ^ 

m .  /Vgo  F 
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(46) 


Adding  (40)  and  (46)  yields  the  total  matrix  element  which  is 
given  by 
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(47) 

and,  then  substituting  (47)  into  (34) ,  one  obtains  an  expres¬ 
sion  for  the  amplitude  dependent  frequency  shift. 
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+ 


1 


(49) 


5.5  Amplitude  Dependent  Frequency  Shift 
a.  Light  Waves 

For  the  circularly  polarized  high  frequency  electro¬ 
magnetic  wave  (o)g  =  a)^e  +  q2c2,  to  >>  a3c j )  /  propagating  par¬ 
allel  to  the  external  magnetic  field,  only  the  last  term  in 
(49)  is  significant  (this  corresponds  to  the  contribution 
from  the  process  shown  in  fig.  26(b))  as  compared  to  the  con¬ 
tribution  from  fig.  24.  Therefore,  expression  (48)  becomes 


Aoo 
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q  J(q)  £ 

8n2 VF2  j  03 3 

o  q  q 


d 3 vf .  (v) 


m .  (03  -qv  ) 

d  q  z 


(50) 


If  the  contribution  to  the  frequency  shift  from  higher  order 
terms  has  to  be  considered,  one  will  have  to  sum  the  series 


where  the  first  term  is  given  by  (50) ,  and  the  beaded  line 
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refers  to  the  screened  Coulomb  potential.  Therefore,  summing 
the  series,  one  obtains 


Hence,  the  modified  expression  for  the  frequency  shift  in  the 
light  wave  will  be 


Aoj_  = 


2  T  (  \  03  4  . 

q Kq)  v  po 

|-  d  3  vf  ^  (v) 

r 

i  — 
8n2VF2  j  (jo  3 

o  q  J  q 

m  .  ( 03  -qv  ) 

J  q  z 

_  _ 

- 1 


(53) 


Since,  for  a  cold  plasma  f (v)  =  n  6 (v) ,  (53)  yields  (see  ap- 

o 

pendix  B) 
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(54) 


where  E  is  the  electric  field  associated  with  the  light  wave. 
The  same  result  can  be  obtained  for  a  wave  polarized  in  the 

counterclockwise  direction. 

Comparing  (54)  with  expression  (50)  of  Montgomery  and 

Tidman  (1964)  one  obtains 
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(55) 


where  Ago  is  the  frequency  shift  in  the  above  reference  which 

is  due  to  the  linearly  polarized  transverse  wave.  Thus,  it 

is  obvious  that  the  result  obtained  in  this  subsection  is 

nearly  the  same  as  that  of  Montgomery  and  Tidman  except  for 

2 


small  factor  of 


8 


1+ 


q2  c2 
3go2 


q 


a 


b.  Whistlers  and  Alfven  Modes 

It  has  been  seen  that  processes  shown  in  figs.  25  and 

26(a)  do  not  give  rise  to  a  frequency  shift  either  for  light 

waves,  whistler  or  Alfven  modes.  However,  the  processes  shown 

in  figs.  24  and  26(b)  do  yield  nonzero  contribution  as  shown 

by  (48)  for  a  cold  as  well  as  hot  plasmas.  Performing  the 

t 

v-integration  in  (48)  for  a  cold  plasma,  one  obtains 
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(56) 


where  the  second  term  is  due  to  ion  contribution  and  the  last 
term  is  due  to  the  process  shown  in  fig.  26(b). 

For  whistlers,  where  go  .  <<  go  <<  go  ,  (56)  yields 

Cl  C  ti 


q 2  I (q) 

8n  VmGO  F 2 

o  q  q 


pe 
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(57) 


Substituting 


I  (q)  = 


into 


(57) 


one  obtains  the  frequency 


VF2 
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shift  for  whistlers 
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e  co 

Pe 

8m  2  co 5 
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(58) 


where  E  is  the  electric  field  associated  with  the  whistler 

2  2 

mode  and  w  =  — -  |  co  |  ,  is  the  linear  frequency. 

q  y2  ce  n  1 

pe 

Similarly,  one  can  obtain  from  (56)  the  frequency 
shift  for  Alven  modes  where  co  <<  co  . 


2  2  2 
q  e  or 
Pe 

16m2  co 5 


CJ 


(59) 


q 


again,  E  is  the  electric  field  in  Alfven  mode,  co  =  qV,  and 

q  A 

V  =  jB2/47Tn  m,  the  Alfven  velocity. 

A  o 

The  frequency  shift  for  electromagnetic  waves  propa¬ 
gating  at  an  angle  0  with  respect  to  the  external  magnetic 
field  can  be  obtained  similarly. 
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APPENDIX  A 


The  nonlinear  dielectric  function  given  by  equa- 

NL  ^ 


tion  (21)  of  chapter  IV  can  be  written 
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where  a  and  a  are  given  by  equations  (45)  and  (46)  of  the 
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same  chapter 
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and  A.  and  Y  ,  Y  are  given  by  (37)  and  (44)  respectively. 

3  1  2 

Since  a  gives  negligible  contribution  as  it  has  been 
i 

seen  from  (48)  and  (49) ,  there  is  no  necessity  for  finding 
—  (e  a  ).  Equation  (A. 3)  yields  after  performing  the  velocity 

o  CO  L  i 
integrations 
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De 


and 


e  (q-q  '  ,  03+03+ ,  )  =  1+ 


1 


03  ‘ 


PI 


|qiq'|2X 


•  12 -v  2 

De 


(  03+03+  ,  )  2 

q  ' 


into  (A. 4),  one  obtains 


1 


-fr- 


Ct  =  — 

2  „2i2 


I 


I  (q'  ) 


q  ASe£L 


n  V  -  ^  „  2 

o  + .  ,+  2mv 
q'^q  e 


(  03+  03  +  ,  )  2 

q 


(o3+03+,  )  2-  |q+q'  |  2c2 


1+  qq' (q+q' ) 

0303^,  ( 03  +  03+ ,  ) 


(g+q 1 ) 

03  +  .  03  +  03+  , 

q  q 


a  +  _ai  + 

0) 


'i  'i  2 


+ 


(  03  -  03^,  )  2 


(  03  -  03 +  ,  )  2-  |q-q'  I  2C2 

q  & 


1+ 


qq'  (q-q1 ) 

0)0)-*-  ,  (  03-03-*-  ,  ) 

q  q 


q  +  ql_  +  q-q' 


q 


31 

2" 

J] 

Multiply  (A. 6)  by  eT  and  differentiate  with  respect 

then  substitute  03  =  qc  and  03+,  =  q'c  to  get 

s  y  ^ 


(A. 4) 


(A. 5) 


(A. 6) 

tO  03  , 
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(q+q* ) 3  +  (q-q1 ) 3 

q2q,2c  04  q2q,2c  04 
So  s  o_ 

(A. 7) 

Performing  the  q ' -integration ,  using  the  spectrum  distribu¬ 
tion  (47) ,  one  obtains 


( e_  a  )  = 
a  GJ  L  2 


16 


q' 


X 2  n  v 
De  o 


I 

q'^q 


Kg') 


X 


mv' 


9 


(e  a  ) 

li  2 


1  1  2 56tt 2  I  (t) 

qCs  q 2 X  2  3N  0 3  mv2 
^  De  D0  e 


(A. 8) 


Differentiate  (A.l)  with  respect  to  co  and  retaining  only 
terms  of  order  I (t) ,  to  get 


^q  '  W  9m  ^Nl/  ^  0)=0)->- 

^  q 


=  0) 


ffL, 

3u)  1  co=(jO^- 

Jj 


go->  ( eT  a  ) 
q  3  oa  L  2 


(A. 9) 


q 


2  oo 2  . 

pi  i 

l 

2  56tt  2 

I  (t) 

r  T  > 
e 

3/2  n 

3 

OO-^ 

-  q 

(qc  )q2X2 
^  s  De 

3n  e3 

D  o 

mv2 

e 

T  . 

^  i; 

Equation  (50)  yields  for  e  =  0. 


Lx)2  . 

PJ- 

2 

OV 

q 


q 2  x 


De 


256tt  2 

I  (t) 

r  T  > 

e 

3/2-7 

1  9N  0 

2 

T. 

L  D  o 

mv 

e 

^  l J 

.. 

(A. 10) 


Substituting  (A. 10)  into  (A . 9 ) ,  one  obtains 


F-* 

q 


q 


2  X 2  (qc  )  ^ 


{( 


i- 


De  s 


256  'j  3/2  128  \ 

—  a)  +  —  a 

30  2  J 


1 


(A. 11) 


and 


a  =  (tt2/0  N  )  ( I  ( t) /mv2  )  (T  /T .  )  3/ 

q  L)  G  v-  j_ 


3/2 


(A. 12) 
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APPENDIX  B 


The  common  ratio  of  the  infinite  series  given  by  equa¬ 
tion  (51)  of  chapter  V  can  be  written  as 


Second  term  in  the  series 
first  term 


(B.l) 


where  the  first  term  is  given  by 


M  = 

and 

Second  term  = 


oo4  q2h2 
pe 


f  (v)  d  3  v 


16n2VF2G02m  ■*  (go  -qv  )  2 
o  q  q  q  z 


(B  .  2 ) 


clu 

^  2  TT  1  > 


Tre2h  }  2  (-i) 


4TTe' 


mVco  F 

q  q- 


h3  4q 2  Ve  ( 2  9 , 2go  ) 

^  q 


y 

f  i2 

HF 

l 

• 

k  ,  k  ,n 
-  x  z 

rE-E  ^  rE+2  -E.  ,  .  ^ 

(  k,n] l  q  k+2q,nj 

oo6  h2q2 

pe  _ 1 _ 

64mn3Va)2F2  e(2q,2oo  ) 

o  q  q  4  q 


d  3  vf  ( v) 

( 00  -qv  ) 
q  z 


( B  .  3  ) 


Substituting  (B.2)  and  (B.3)  into  (B.l),  one  obtains 


GO 


pe 


4n  e  (2q , 2go  ) 
o  4 


d 3 vf (v) 

(go  -qv  ) 
q  ^  z 


(B .  4  ) 


This  yields  for  a  cold  plasma 


GO' 


= 

4co2 

q 


e  (  2q  ,  2go  ) 

q 


(B.  5) 


. 
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and 

o 

go"  , 

pe  1 

4go2  e  (  2q  ,  2go  ) 

q  n  »■  q 

where  the  dielectric  function  for  light  waves 


£(20,200  ) 

q 


4oo  2 

q 


Substituting  (3.7)  into  (B.6)  yields 


3  1+q2c2//3a)q 

2  1+q  2  c  2  /go  2 

q 


(B.6) 


is 

(B .  7 ) 


( B  .  8  ) 


B30022 


